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ABSTRACT 


The aging E-2C fleet is expected to be retired by the year 2015. In order to 
orovide Airborne Early Warning (AEW) for the battle group during the 
transitional years and beyond, the design of a replacement aircraft must begin 
soon. In order to conform with present day economic realities, one possible 
configuration is a new airframe using the radar system and rotodome which 
currently operates on the E-2C. Other likely requirements for a new AEW 
aircraft includes a high-speed dash (M=0.7-0.85) capability, an extended 
mission time (up to 7.5 hours), turbofan engines, and an aircrew ejection 
system. 

The results of this design effort includes an investigation of a possible 
configuration and the aerodynamics involved. Performance and Stability & 
Control characteristics are also discussed briefly. Finally, a qualitative analysis 


of the use of the E-2C's radar system on a new airframe will be presented. 
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l INTRODUCTION 


The purpose of this thesis is to provide an initial conceptual design for a 
carrier-based Airborne Early Warning (AEW) aircraft that would replace the 
E-2C. The AEW aircraft design is in response to a Proposed Request For 
Proposal (Proposed RFP), which is based on the perceived need to replace the 
E-2C. The Proposed RFP was prepared by C.F. Newberry after informal 
discussions with several individuals including students, Naval Air Systems 
Command (NAVAIRSYSCOM) staff, and other members of the E-2C community. 
It is not an official document, but rather a general guideline for an AEW design. 
The Proposed HFP is included as Appendix A. This chapter will provide some 
introductory material necessary to understanding the issues involved in 
designing any generic AEW aircraft. A description of a generic AEW mission 
profile will be discussed. Additionally, a brief description of the method of 


design will be presented. 


A. BACKGROUND 
1. Proposed Request For Proposal 
With an increasingly aging E-2C fleet, the Navy has recently 
recognized the need for a replacement AEW aircraft. In accordance with 
present economic realities, the first objective is to provide a capable platform 


that is cost effective. A “low risk airframe configuration” is most desired. A low 


risk detection system is also desired. In order to satisfy the above objectives, a 
Proposed RFP requirement is to include the existing 24-foot rotodome currently 
being used on the E-2C in the new design. 

In order to detect high-speed adversary aircraft as far from the battle 
group as possible, and to quickly replace an aircraft with an inoperative 
detection system, there is a requirement that a new AEW platform possess a 
high speed dash (M=0.70-0.85) capability. The aircraft must also possess 
excellent loiter characteristics in order to provide long periods of detection for 
the battle group. A total unrefueled mission cycle time of 5.75 hours is required. 
Additionally, an in-flight refueling capability is required to extend mission cycle 
time. 

The new AEW aircraft is required to provide direct self defense. It is 
expected that two AIM-7 Sparrow-sized missiles would be mounted on wing 
Stations. Additionally, it is required that the aircraft possess chaff and flare 
launchers. Also, there is a requirement for a crew ejection escape system. 

Carrier Suitability requirements include total compatibility with all 
CVN-68 (Nimitz class) carriers and subsequent, and a maximum takeoff weight 
of 60,000 Ibs. Also, in an effort to remove the hazards of spinning propellers on 
the flignt deck, a turbofan propulsion system is required. Table 1 outlines the 
significant Proposed RFP requirements for the AEW aircraft. 

2. AEW Mission Profile 
The Proposed RFP specified some general mission requirements the 


AEW aircraft must be able to accomplish. Also included is standard information 


on essential mission parameters such as start, taxi, fuel reserves, etc. These 
requirements were used along with a baseline knowledge of the AEW mission 
to generate the mission profile shown in Figure 1. Mission parameters are 


summarized in Table 2. 


ПОЕНЕ НО БЕО АБР БЕО БНЕМЕМТ5 


PROPOSED RFP РІС REQUIREMENT 
High Speed Dash Mach = 0.70-0.85 
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Figure 1. AEW Mission Profile 


It should be noted that some of the performance parameters presented 
in the Mach number, Distance, and Time columns in Table 2, are approximated 
based on historical trends and past experience. A more detail estimation of 


performance is provided in Chapter V. 


TABLE 2. MISSION PARAMETERS 


PHASE ALTITUDE | DIS- TIME | TOTAL | POWER 
NO | (ЕТ) | TANCE TIME 
LM E m че е 
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Also note that by choosing a specific Mach number for the high speed 














dash phase, the first design decision was made. The Mach number range 
given in the Proposed RFP was too broad. The upper end of the Mach number 
range seemed a little too high (M=0.85), particularly from the standpoint of drag 
divergence. On the other hand, the lower end of the range (M=0.70) seemed a 
little too low from the standpoint of design technology. It was decided that a 
mid-range Mach number (M=0.78) was the maximum realistic speed to which 


this AEW aircraft could be designed. 


B. DESIGN STRATEGY 

As previously mentioned, the primary purpose of this research was to 
provide a first iteration on a conceptual design only. As such, the areas of 
research are directly proportional to the areas of emphasis given in the 


Proposed RFP. The focus of this research will be on the aircraft configuration 


and the resulting aerodynamics. Performance and Stability & Control will also 
be discussed briefly. Some of the topics addressed in preliminary design books 
such as References (1) and (2) are outside the scope of this research. Such 
topics include propulsion, structures, and cost analysis. A more complete 
design effort is possible only after an entire design team is assembled. 

The primary objective during the design process was to remain focused on 
what the customer (NAVAIRSYSCOM) might desire in a AEW aircraft. This 
design approach, known as Quality Function Deployment (QFD), seems 
obvious but is a new concept to most design teams. QFD will be discussed in 
detail in Chapter Il. 

In order to avoid “reinventing the wheel” and to keep costs down, 
characteristics of proven aircraft with similar missions (i.e., E-2C, S-3A, EA-6B) 
were evaluated, and integrated into this AEW aircraft design. The overall 
philosophy was to keep the AEW aircraft design as simple, and as conventional 
as possible. Design techniques and equations were used in accordance with 
conventional design books such as References (1) and (2). Also, computer 
programs such as MATLAB and EXCEL were used as much as possible to 
rapidly complete future iterations. The programs are included as appendices. 
The equations in each computer program are referenced with the appropriate 


book and equation number, in order to assist any follow-on work to this thesis. 


Il. PRE-DESIGN ANALYSIS 


It is widely understood that the further along a product is in its design 
process, the less design freedom the engineer enjoys. Therefore before any 
design process begins, it iS imperative that the customer’s desires and 
parameter constraints be thoroughly analyzed. This chapter will examine the 


specifics of QFD, and the constraints placed on the AEW aircraft. 


A. QUALITY FUNCTION DEPLOYMENT (QFD) 

Because of the present realities of fierce global competition, major 
companies throughout the world are searching for creative ways to produce 
high quality products at competitive prices. For governments on tight budgets, 
the commitment to high quality and low cost has also become increasingly 
important. The results of these realities have been numerous quality-based 
management, engineering, and design philosophies. Some of these 
philosophies include Deming's Total Quality Management (TQM), Taguchi's 
Parameter Design Method, and Mitsubishi's Quality Function Deployment 
(QFD). It has been these kinds of quality-oriented philosophies that have made 
Japanese industries so successful. Because these strategies are 
complementary, the more general term of QFD will be used for the purpose of 


this discussion. 


As noted in Reference (3), it is extremely difficult (and costly) to implement 
quality into a product that has already been designed. Therefore in order to 
design a quality product, it is imperative that before a preliminary design 
process begins, sufficient time must be spent on the issue of product quality. 
From the standpoint of QFD, the answer to the question “What is Quality?” is 
simple--quality is providing what the customer wants! Reference (4) provides a 
more formal definition--“Quality is the loss a product causes to society after 
being shipped, other than any losses caused by its intrinsic functions”. The 
purpose of QFD is to investigate what the customer wants in detail, and then 
translate those desires into engineering and design decisions. 

The result of implementing QFD speaks for itself. As Reference (5) points 
out, Toyota Auto Body reduced costs by 61% after implementing QFD. 
Reference (6) notes that an unspecified Japanese automaker with QFD takes 
32 months from first design to finish a car, while it takes 60 months for a U.S. 
automaker without QFD! These results were accomplished because of a 
commitment to begin the design process only after extensive customer research 
was completed. Once the design process was underway, the need for design 
changes became almost non-existent, because the customer's desires were 
already known. Figure 2 is reproduced from Reference (5) and graphically 
illustrates the difference in the design philosophies between two automobile 
companies. The lesson to be learned is clear--if more time and money are 
spent investigating customer desires before the design process begins, more 


time and money will be saved in the long run, and product quality will be higher. 
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Figure 2. Results of QFD [Ref. 5] 

In terms of an AEW aircraft design, a preliminary QFD analysis was 
performed based on the customer's (NAVAIRSYSCOM's) perceived desires 
expressed in the Proposed RFP. These desires, commonly referred to as 
Customer Attributes (CAs), were then numerically prioritized in accordance with 
the relative importance E them in the Proposed HFP. Based on the 
customer attributes and their relative importance, a House Of Quality (HOQ) was 
constructed. The HOQ is a matrix-type figure that puts customer attributes into a 
format that is usable by both engineering and management. The HOQ is shown 
in Figure 3. 

Several items should be mentioned in the construction and use of the 


HOQ. As was previously mentioned, CAs were ranked according to the relative 


importance given them in the Proposed RFP. The Relative Importance (RI) is an 
integral part of the HOQ because it is a constant reminder to both management 
and engineering of their priorities. The RI is a major tool for making design 


decisions. 
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Figure 3. House of Quality 
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Note that Figure 3 shows CAs vs. Engineering Characteristics (ECs). The 
CAs can be considered the “what” portion of the HOQ while the ECs can be 
thought of as the “how” portion. This is because the CAs communicate what 
needs to be accomplished while the ECs tell us how they can be 
accomplished. Reference (5) points out that, “Engineering Characteristics 
should describe the product in measurable terms and should directly affect 
customer perceptions”. Thrust-to-Weight ratio (T/W) for example, is clearly 
measurable and it will directly affect how the customer perceives the product in 
terms of its performance characteristics. Also note that shown with each EC is a 
plus or minus sign. This communicates to the engineer what should ideally be 
accomplished with a particular EC. For example, the Weight EC is followed by a 
minus sign because the objective is to keep weight as low as practical. 

The central matrix portion of Figure 3 is the primary vehicle in which CAs 
and ECs communicate. As Reference (5) notes, it is in this central matrix that 
ECs that affect particular CAs are identified, and relationships between them 
are established. For example, there is a positive relationship between low 
Weight (EC) and maximum Endurance loiter (CA). In other words, all other 
things being constant, the lower the weight the longer the loiter time. Once this 
matrix is completed, the engineer will have a better idea of how to proceed in 
terms of the design process. 

Another significant part of the HOQ is the characteristic roof. The roof is 
used to establish relationships between various ECs. For example, there is a 


negative relationship between low weight and higher Fuel Volume. Like the 
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central matrix, the completed roof helps the engineer make the necessary 
decisions in the design process, by balancing these relationships. 

The HOQ shown in Figure 3 is only the first in a series of four or more 
HOQs that can be used to communicate the customer's desires through to the 
actual manufacturing process. Figure 4 is reproduced from Reference (5) and 
shows an example of how these HOQs might be related and how CAs trigger a 
series of decisions made through to manufacturing. Note that the "how" portion 
of each HOQ becomes the "what" portion of the next HOQ. The subsequent 
HOQs in the series would necessarily be generated after future iterations in the 
design process. It is difficult for example, to examine the characteristics of 


specific parts while still in the conceptual phase. 





Figure 4. Linked HOQs [Ref. 5] 
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It should be emphasized that the HOQ shown in Figure 3 is preliminary. It 
is based on the preliminary requirements given in the Proposed RFP,and is 
primarily used for setting design priorities. Before the AEW aircraft design goes 
beyond the conceptual phase, detailed marketing research should be 
conducted to investigate what the customer wants. The research should 
include a survey of all the customers including NAVAIRSYSCOM, aircrew, and 
maintenance personnel. The research should be a study of likes and dislikes of 
even the smallest details of an AEW aircraft. For example, questions on the 
operation of the external door, or the location of a parking brake, etc., should be 
included when questioning customers. This research would then generate 
many series of HOQs. 

The QFD strategy cannot be overemphasized in the aircraft design 
process. Although the process may seem time consuming and wasteful at first, 
a properly implemented QFD program will result in enormous long run benefits 
to both the aircraft company and the customer. Within the scope of this 
research, only aircraft companies with fully implemented QFD programs should 


be considered for development of the AEW aircraft. 


B. CONSTRAINT ANALYSIS 

Before the actual design process can begin, it is necessary to evaluate two 
of the aircraft's characteristics. These characteristics are T/W and Wing 
Loading (W/S). A series of performance equations may be derived in which 


T/W is expressed as a function of W/S. These equations are derived in 
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Reference (7). Equation constants are obtained from performance 
characteristics provided in the Proposed RFP. For a range of W/S, a range of 
T/W may be generated for each equation. The equations are then graphed ona 
single constraint plot. The plot graphically depicts a solution space. Any T/W- 
W/S combination may be selected within that space. Obviously, some T/W-W/S 
combinations will be better than others. For example, suppose a constraint 
analysis On an aircraft reveals that lowest T/W in the solution space is 0.25. 
This means the aircraft can perform the required mission at a T/W = 0.25. It 
would be illogical to choose a T/W = 0.50 even though it is also within the 
solution space. It should be noted that although the constraint plot is primarily a 
pre-design tool, it may be used throughout the design process. As more 
knowledge of the design is known, more exact iterations of the constraint plot 
may be generated. It should also be pointed out that the constraint analysis 
need not be limited to performance equations only. For example, if a valid 
expression for maintainability in terms of T/W and W/S is found, it should also 
be included as part of the constraint analysis. 

In order to keep future iterations simple, a computer program was written in 
MATLAB, based on the performance equations derived in Reference (7). The 
complete program is included as Appendix B. All equations in Reference (7) 
applicable to the AEW mission were used with the exception of takeoff and 
landing performance. Expressions presented in Reference (1) were used for 
takeoff and landing performance because of their simplicity and their more 


conservative results. Performance equation constants were obtained from 
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performance characteristics provided in the Proposed RFP and from a baseline 


knowledge of the AEW mission. The results of the AEW constraint analysis is 


shown in Figure 5. 
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Figure 5. AEW Constraint Analysis 


The solution space is the outlined upper center portion of the graph. Note 


the relatively flat bottom of the solution space. This flat bottom is most fortuitous 


because it allows a certain degree of design freedom. For a relatively low 


S 


T/W of 0.46, a W/S anywhere between 55 and 116 Ibs/ft2 can be chosen. 
Because of wing area limitations for carrier operations however, the W/S for an 
aircraft of this size is typically between 70 and 116 Ibs/ft2. 

Also note that the constraint plot includes a maintainability line. The line is 
the result of a equation derived in an unpublished paper by C.F. Newberry. The 
equation is the result of a linear curve fit of data from 25 different aircraft. | 
should be noted that there are limitations in the application of this equation. 
First, none of the aircraft for which data was supplied are Navy aircraft. Navy 
aircraft traditionally have different Mean Man Hours/Flight Hour (MMH/FH) rates 
than other aircraft. Second, a general trend should not be assumed using 25 
very different aircraft. These aircraft ranged from T-38's to 747's. Although the 
validity of the maintainability line may be suspect, it should be investigated in 
greater detail, using a larger database of aircraft similar to the aircraft being 
designed. The current maintainability equation may be used in the constraint 


analysis, but only as long as its impact is integrated in a reasonable fashion. 
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lll. АЕМ CONFIGURATION 


This chapter will discuss the initial conceptual design for the AEW aircraft. 
A description of the aircraft will be provided along with the rationale behind 
various design decisions. An initial weight & balance evaluation will also be 
discussed. Finally, an analysis of the AEW aircraft with various carrier suitability 


requirements will be performed. 


A. AIRCRAFT DESCRIPTION 
1. Introduction 
The purpose of this section is to provide a brief description of the 
external aircraft configuration, and to provide justification for some design 
choices. Not all configuration characteristics of the aircraft will be discussed in 
this section however. Aircraft characteristics directly related to aerodynamics 
will be discussed in Chapter IV. These characteristics include planform 
Selection, airfoil selection, and high lift devices. 
2. General 
The AEW aircraft design is shown in Figure 6. The aircraft is designed 
to hold a crew of four and will be powered by twin turbofan engines. Crew 
seating will be arranged in a dual-tandem configuration. Large cockpit 
windows will allow better visibility for carrier (CV) launch and recovery 


Operations. The rotodome antenna will be supported by the existing rotodome 
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Figure 
6. AEW Aircraft Design 


pylon. Also, in order to satisfy CV requirements, the rotodome retraction system 
that was operational on early E-2's must be used. Twin vertical stabilizers will 
be mid-mounted at either end of the horizontal stabilizer. A total fuel weight 
estimate of 14000 pounds was based on fuel volume calculation procedures set 
forth in Reference (8). It should be noted that this iteration of the aircraft design 
includes no composite materials. Significant aircraft dimensions are presented 
in Table 3. 
3. Specific Component Description 
a. Engines 

Although a detailed study of the propulsion system was outside 
the scope of this design effort, an initial analysis of the required engine 
performance was made. In order to meet the mission requirements of high- 
speed dash and long time loiter, it is clear that a high-bypass turbofan engine 
with a low Thrust Specific Fuel Consumption (TSFC) is required. Assuming an 
initial takeoff weight of approximately 55,000 los. and a T/W = 0.46, the thrust 
per engine requirement is approximately 12,700 lbs. As shown in Reference 
(9), the technology for such an engine already exists. Two operational engines 
with characteristics similar to those required for the AEW aircraft, are presented 
in Table 4. Further design iterations should include an investigation into the 
feasibility of using an upgraded version of the General Electric (GE) TF34-GE- 


400A engine in the AEW aircraft. 
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TABLE 3. AEW AIRCRAFT DIMENSIONS 
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Elevator Area ЗИ 2 
Vertical Tail Area 90 ft2 
Vertical Tail Sweep 26.6 degrees upper, 36.9 degrees 
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TABLE 4. SIMILAR ENGINE CHARACTERISTICS 


Engine | Maker | Type Thrust | TSFC | Pressure} Dimen- | Weight 
(105.) 1 Капо | sions (Ibs.) 

| | | wes 

TF34-GE1 General |AFF з 9,275 07863 EN S2in. x | НЕМА 


FJR-710-|Nat. Aero.| AFF 3 14,330 | 0.340 22 57.lin. x] 2 
/600S 4 Lab 92.5in. 
Tokyo 


Notes: 1- At Maximum Power 
2- S-3A Aircraft 
3- Axial Flow Fan 
4- NAL/Kawasaki Aircraft 
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The engines should be mounted closely to the wing for two 
reasons. First, exhaust flow through the slotted trailing edge flaps will help 
reattach the airflow over the wing, thereby increasing CLmax. Second, an 
engine mounted closely underneath the wing is further from the ground, and 
therefore less likely to ingest foreign objects. This would result in fewer engine 
replacements and lower life cycle costs. 

b. Vertical Tail 

As previously mentioned, the empennage will include two vertical 
Stabilizers. The maximum height of the vertical stabilizers were modeled after 
the E-2C in an effort to keep the tails from interfering with the look-down 
capability of the rotodome antenna. Each vertical stabilizer will include a rudder 
control surface. It should be noted that if future iterations mandate higher 
vertical tails, maximum use of composites will be necessary to avoid antenna 
interference. 

c. Aircraft Entry 

Aircraft ingress will be accomplished through a single door in the 
fuselage. A walkway will allow movement between the door and the cockpits. 
The major advantage of this configuration is flexibility. The walkway will allow 
the crew to move freely throughout the aircraft to troubleshoot avionics systems, 
Switch seats, etc. Consideration may be given to a canopy system similar to that 
currently operating in the EA-6B. The canopy arrangement was initially ruled 
out in this study due to potential engineering difficulty, increased life cycle costs, 


and lack of flexibility. 
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d. Wing Fold System 

The first wing fold will be at 15 feet from the aircraft centerline. 
This will result in a maximum wing fold span of 30 feet. This wing fold span is 
within the maximum requirement of 35 feet and will allow easy storage of 
aircraft on the flight deck. The wings are intended to fold vertically up. At the 
completion of this vertical fold, the wing tip will physically interfere with the 
rotodome antenna. Therefore a second wing fold at 30 feet from the centerline 
is required. Dashed lines denote the wing fold breaks in Figure 6. The 
horizontal wing fold system which currently operates on the E-2C was ruled out 
for two reasons. First, horizontally folded wings create a large sail area. When 
the aircraft taxis perpendicular to the wind on the carrier deck, it tends to get 
blown, resulting in lose of control. Second, it is clear from the geometry of this 
AEW design that the wingtip of a horizontally-folded wing would not reach a 
wing support on the horizonal tail tip. 

e. Armament 

The aircraft is designed to accommodate one wing station on each 
wing at approximately 14 feet from the centerline. Each wing station should be 
capable of carrying an air-to-air missile of 500 pounds. Although use of the 
AIM-7 Sparrow missile was alluded to in the Proposed RFP, this is not 
recommended. Use of the AIM-7 would require the aircraft to possess a high- 
energy, target illumination capability. The new generation of "fire-and-forget" 


air-to-air missiles such as AMRAAM and Have-Dash are much more suitable for 
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the AEW aircraft. No target illumination is required for these missiles. Updated 
target information is provided via data link. 
f. Landing Gear 

A landing gear analysis was performed based on procedures set 
forth in Reference (2). The aircraft will use a standard tricycle system. 
Longitudinal placement of the main gear was determined by an estimated 
center of gravity location. Lateral placement of the main gear was determined 
by a Maximum overturn angle requirement of 54 degrees. The wheelbase will 
be 26 feet long and the main wheel width will be 20 feet. The nose gear will 
have a dual-wheel configuration. The nose gear will retract aft into the 
fuselage. Each of the main landing gear will be a single-wheel configuration 
and will also retract aft into the fuselage. Approximate tire dimensions are 25 in. 
x 7 in. (diameter x width) for the nose and 45 in. x 17 in. for the main. These 
dimensions are approximately 25% greater than the statistical equation 
proposed by Reference (2). This dimensional increase is to account for the 
harsh landing environment of the aircraft carrier. The 25% dimension increase 
corresponds well with the tire sizes of current carrier aircraft. 

g. Escape System 

The Proposed RFP requires the installation of an all-crew ejection 
system in the AEW aircraft. This requirement has resulted in many difficulties in 
the design of the escape system. These difficulties are obviously the result of 
the rotodome. An approximate trajectory of the aircrew on ejection is shown in 


Figure 7 for three flight conditions. An ejection trajectory computer program was 
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written in MATLAB and is included as Appendix C. The parabolic 
approximation is based on an ejection analysis presented in Reference (10). 
The identical pair of trajectories represent the front seat and back seat ejections. 
The diamond figure represents the location of the rotodome antenna. 

It is obvious from the Figure 7 that the ejection system will result in 
aircrew impact with the rotodome. A bottom or sideways ejection would require 
development of a new ejection system, and obviously could not provide a 0/0 
ejection capability. After an examination of various aircrew and rotodome 
placements, it became apparent that with today's technology, there are no safe 
ejection alternatives with the rotodome installed. 

Ejection of the rotodome prior to crew ejection also has significant 
problems. [he rotodome antenna alone (not including the supporting pylon 
and shaft) weighs 2350 pounds. In order to get the crew out of the aircraft 
quickly, the rotodome would have to be ejected with a typical acceleration of 
approximately 12g's. This would require a series of rockets that would have to 
generate a combined force of over 28000 pounds. These rockets would most 
likely have to be very large in order to provide such a force. It is unlikely that the 
rockets would fit into a supporting pylon that is only approximately one foot 
wide. 

Additionally, it is obvious that the rockets would have to be directly 
attached to the rotodome. This means they would rotate with the rotodome. 
This means there would be no way to direct the trajectory of the rotodome, 


because it must be ejectable at any time during the rotation. Therefore, the 
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rockets would have to be of equal propulsive force. During certain flight 


conditions, including a 0/0 ejection, the crew would still be in danger of ejecting 


into the rotodome. 
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Figure 7. Aircrew Ejection Trajectory 


Ejecting the entire rotodome structure would eliminate the 
controlled trajectory problem, but would generate other problems. Now the 
rockets would have to generate a combined force of over 38000 pounds. The 


rockets under the forward supports would most likely ignite the fuel in the fuel 
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cells directly below. The resulting explosion would jeopardize the lives of the 
aircrew during ejection. 

Two final points are worth mentioning. First, the new technology 
and the resulting developmental costs of ejecting a rotodome will likely be 
enormous. Second, any further investigation into rotodome ejection should 
necessarily include an examination of how the pitching moments about the 


center of gravity are affected . 


В. WEIGHTS, CENTER OF GRAVITY, AND MOMENTS OF 
INERTIA 
1. Weights 
An evaluation of the AEW aircraft weight was performed using the 
individual component equations given in References (1) and (8). A computer 
program was written on MATLAB using the applicable equations. Many of the 
equations represented individual weight components as a function of takeoff 
weight. Since the determination of the takeoff weight was the ultimate objective, 
the program uses a secant method iteration procedure to find the takeoff weight. 
The weight program is included as Appendix D. In order to assure the accuracy 
of the program, a weight analysis on the E-2C was performed. It was found that 
the program prediction came within 300 pounds of the actual E-2C weight. The 
program was then used to analyze the weight of the AEW aircraft. The 
predicted weight was found to be approximately 53000 pounds which is 


comparable to the E-2C weight and well within the maximum requirement of 
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60000 pounds. The aircraft possesses a 7000 pound weight growth potential 
for future avionics upgrades. 
2. Center of Gravity and Moment of Inertia 

Component weights calculated from the weight program were used to 
approximate the aircrafts Center of Gravity (CG) and Moment of Inertia. 
Component CG locations were approximated based on procedures set forth in 
References (1), (2), and (8). Component Moment of Inertia values were 
calculated in accordance with procedures set forth in References (2). The 
component characteristics were used to calculate aircraft CG and Moment of 
Inertia values. All calculations were performed on a computer program written 
on EXCEL. The computer program was acquired from Reference (11). The 
computer program and the results of this program are included as Appendix E. 
An initial approximate CG location is 32.4 feet aft from 5 forward of the nose 
(approximately 48.6% MAC), and 10.9 feet up from 5 feet below the fuselage. 
More detailed CG and Moment of Inertia calculations will obviously be 


necessary with future iterations of the design. 


C. CARRIER SUITABILITY REQUIREMENTS 
Carrier suitability dimensional requirements and the significant AEW 


aircraft dimensions are shown in Table 5. 
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TABLE 5. CARRIER SUITABILITY DIMENSIONAL COMPARISON 


DIMENSION REQUIREMENT AEW AIRCRAFT 
Max. Gross Weight 60000 Ibs. 53000 Ibs. 
Max. Wing Span 








Max. Height 15.5 fh 18.5 ft. (rotodome 
retracted) 
Max. Main Gear Width 






Min. Tipback Angle 20 deg. 
Max. Tipover Angle 52.5 deg. 
Elevator Size Restriction 52 X 85 ft. 55 X 9990 
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IV. AERODYNAMICS 


In order to get maximum effectiveness from an airframe and its propulsion 
system, a thorough examination of the aircraft's aerodynamic characteristics 
during the design process is mandatory. This chapter will examine the design 
decisions involved in selecting the AEW aircrafts airfoil and wing planform. 
Additionally, the aircraft's lift curve slope and high lift devices will be discussed. 


Finally, an analysis of the aircraft's drag characteristics will be presented. 


A. AIRFOIL SELECTION 

Because of the Proposed RFP requirements, the AEW aircraft will be 
expected to operate under a variety of flight conditions. It must be able to cruise 
at high subsonic speeds, loiter for long periods of time, and possess carrier- 
suitable, slow flight characteristics. In order to meet these requirements, the 
wing's airfoil must possess several seemingly contradictory characteristics. 
The airfoil should have a relatively high thickness ratio in order to increase 
Clmax, increase benefit from high lift devices, decrease weight, and increase 
wing fuel storage capacity. If the wing is too thick however, the drag divergent 
Mach number (Maa) will be too low to satisfy the high speed dash requirement. 
An increase in Mgg could be accomplished through an increase in wing sweep, 
but this generates additional problems which will be discussed in the next 


section. The airfoil must also have a high Clmax for the loiter and landing 
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phases of flight. Most high speed airfoils however, are not known for their high 
Clmax values. Finally, the airfoil’s thickness distribution should be investigated 
in terms of its skin friction drag characteristics. As Reference (12) notes, a 
maximum thickness that is close to the trailing edge results in a more favorable 
pressure gradient on the forward portion of the airfoil. This helps create more 
laminar flow which results in reduced skin friction drag. It should be noted 
however, that an aft maximum thickness can cause poor pressure recovery 
characteristics at high angles-of-attack. 

Based on the above requirements, it became clear that a supercritical 
airfoil was necessary. A Supercritical airfoil is characterized by a relatively flat 
upper surface, and a maximum thickness located near the trailing edge. It also 
has a relatively blunt leading edge, and it is cambered at the aft portion of the 
airfoil. Reference (13) notes that for a given thickness ratio, the supercritical 
airfoil has a higher Mgg than conventional airfoils. This allows a thicker wing 
and less wing sweep. Additionally, the Supercritical airfoil has a much higher 
Clmax than a comparable conventional airfoil. Finally, the thickness distribution 
and the trailing edge upper and lower surface tangency results in a more 
favorable pressure gradient. The aft maximum thickness of the supercritical 
airfoil does not result in pressure recovery problems, because the camber is 
accomplished primarily by the lower surface. This allows the upper surface to 
remain relatively flat. 

It should be pointed out that use of a supercritical airfoil will not be without 


its difficulties. First, the very thin trailing edge could prove to be a structural and 
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manufacturing problem. Second, although the original supercritical airfoil was 
designed in 1965, development and testing of an entire family of supercritical 
airfoils has been relatively recent. Because Supercritical airfoils are relatively 
new technology, development costs may be high. Finally, the aft camber of the 
airfoil will result in large negative pitching moments. Despite the potential 
difficulties however, the supercritical airfoil shows the most promise in terms of 
Satisfying the requirements of the Proposed RFP. 

Initially it was hoped that an airfoil with a thickness ratio of 0.14 could be 
used for on the aircraft. Even with some compromise in the wing sweep, it soon 
became evident that a lower thickness ratio would be necessary in order to 
reach an acceptable Mgg. Experimental data presented in Reference (14) 
shows that at a thickness ratio of 0.12 and a design Cl of 0.7, the airfoil Mag is 
approximately 0.76. A moderate wing sweep should permit reasonably low 
drag characteristics at the design cruise Mach number of 0.78. 

After an evaluation of the family of NASA supercritical airfoils, it became 
clear that the best airfoil for the required mission was the NASA SC(2)-0712. 
This airfoil is shown in Figure 8. The airfoil's coordinates are reproduced from 
Reference (14), and is included as Appendix F. An explanation of the NASA 


supercritical airfoil designation system is presented below. 
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Supercritical phase 2. Design lift Thickness 
There are currently 3 coefficient Ratio 
phases of airfoil designs. (tenths) (hundredths) 
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One of the biggest difficulties in selecting an airfoil was in obtaining the 
specific airfoil characteristics. Because of the relatively new technology, there is 
no compiled source of information for supercritical airfoils (such as Reference 
(15) for conventional airfoils). The three sources that provided most of the 
information on the airfoil were References (14), (16) and (17) . Airfoil 


characteristics are presented in Table 6. 
0.3 
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Figure 8. NASA SC(2)-0712 Airfoil 


TABLE 6. NASA SC(2)-0712 CHARACTERISTICS 
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B. PLANFORM DESIGN 

Given the target cruise Mach number of 0.78 and the relatively thick airfoil, 
it was clear a planform with significant wing sweep would be required. Too 
much wing sweep however, generated numerous problems including a 
decrease in Clmax and CL., increased wing weight and decreased wing fuel 
volume. Selection of the previously mentioned airfoil was made only after it 
was determined that a relatively high Mag could be attained with a modest wing 
sweep. 

Figures 9 and 10 show the results of trade studies conducted to graphically 
illustrate the parameters involved in planform design and airfoil selection. 
Figure 9 shows Mgg as a function of thickness ratio with varying sweep. Figure 
10 shows how thickness ratio and wing sweep affect wing weight. The results 
of these parametric studies were used to select the optimum planform design 
and airfoil thickness. With an airfoil thickness ratio of 0.12, a leading edge wing 
sweep of 21 degrees is the optimum choice considering all the parameters 
involved. This results in a wing Mgg of 0.81. 

With the leading edge wing sweep selected, the focus of attention was then 
directed to the trailing edge sweep. A trailing edge sweep of 6.5 degrees was 
selected for a first iteration. The relatively small sweep will insure efficient use 
of flaps and aileron control surfaces. The flatter trailing edge sweep also allows 
an increase in wing area and wing fuel volume. With a wingtip chord length of 
four feet selected as a first iteration, and the above planform characteristics, a 


wing area of 639 ft2 was calculated. 
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Figure 9. Wing Mgg With Varying Wing Geometry 
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Figure 10. Wing Weight With Varying Geometry 
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Another consideration in the planform design was aspect ratio. It was clear 
that in order to satisfy aggressive loiter requirements, a high aspect ratio would 
be necessary. For a given wing area, this would mean a larger wing span. Too 
large a wing Span causes two problems however. First, it would result in line-up 
difficulties during carrier landings. Second, the large wing span would result in 
Signal interference with the rotodome antenna, degrading radar performance. 
The selected wing span of 72 feet results in a aspect ratio of 8.11. The resulting 


maximum L/D ratio is 16. 


C. LIFT CURVE SLOPE 

With the selection of the wing planform design, a calculation of the wing’s 
lift curve slope was then possible. Calculations were done in accordance with 
the procedures set forth in References (1), (2) and (18). The lift curve slopes for 


three flap settings are shown in Figure 11. 


D. HIGH LIFT DEVICES 
In order to make landing speeds slow enough to meet the Proposed RFP 


carrier Suitability requirements, а СЇ тах of approximately 3.0 is required. To 
accomplish this, double slotted flaps are necessary. In accordance with the 
procedures set forth in Reference (2), ACLma, and Aeg values were calculated. 
А тахїтит А СЇ лах was calculated to be 0.98. 

Two design characteristics that will help increase CLma, with the flaps 


down should be mentioned. First, engines should be situated on the wing so 
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that engine exhaust will flow through the slotted flaps. Second, use of a aileron 


droop system with the flaps will help increase the CL a, of the entire wing. 
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Figure 11. AEW Lift Curve Slope 


E. PARASITIC DRAG CALCULATION 


Parasitic drag (CDo) calculations were performed in accordance with 


procedures set forth in Reference (18). A CDs computer program was written in 


MATLAB and is presented in Appendix G. A CD, of approximately 0.0205 was 
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computed by the program. This CD, value will be used to calculate a drag polar 


for the AE W Aircraft. 


F. DRAG POLAR 

The AEW drag polar was computed assuming CD as a parabolic function 
of CL. A first iteration efficiency factor of 0.8 was assumed. Also, the previously 
determined aspect ratio of 8.11 and CD, of 0.0205 were used in the equation. 
A drag polar for the AEW aircraft in the clean configuration is shown in Figure 


2; 
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Figure 12. AEW Drag Polar 
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V. PERFORMANCE 


This chapter will present the results of a preliminary performance analysis 
conducted for the AEW aircraft. This analysis was primarily performed using a 
computer program written in MATLAB. The program is presented in Appendix 
H, and also includes some aerodynamic calculations such as Coefficient of 
Drag (Cp) and Lift-to-Drag ratio (L/D). A Takeoff and Landing computer 
program is also included in Appendix H. Performance calculations were done in 
accordance with References (1) and (19). The equations in the programs are 
denoted with the equation number from the appropriate Reference. For all 
performance characteristics, it has been assumed standard day unless 
otherwise noted. Additionally, all results were generated for the clean 
configuration, with the obvious exceptions being the takeoff and landing phases 


of flight. 


A. Takeoff and Landing 

Because of the angle between the aft landing gear, the vertical stabilizers 
and the ground (see Figure 6), it is necessary to limit aircraft rotation to no more 
than 18 degrees. This angle of rotation is sufficient however, because the 
typical rotation on takeoff is approximately 10 degrees. References (1), (2) and 
(19) provided schematics and distance equations necessary for takeoff and 


landing. Takeoff and landing schematics are shown in Figures 13 and 14, and 
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are reproduced from Reference (1). Takeoff and landing distances are shown 


in Tables 7 and 8. 
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Figure 13. Takeoff Schematic [Ref. 1] 


TABLE 7. TAKEOFF DISTANCES 
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Figure 14. Landing Schematic [Ref. 1] 


TABLE 8. LANDING DISTANCES 


Landing Distances Standard Day 












Hot Day (90°F) 





B. Thrust Required 

The thrust required for the AEW aircraft at three altitudes between sea level 
and 35,000 feet are shown in Figure 15. The calculated thrust required curves 
were used to generate other performance characteristics such as power 


required and rate of climb. 
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Figure 15. AEW Thrust Required 


C. Power Required and Power Available 

AEW Power Required and Power Available Curves at sea level, 15000 ft, 
and 35000 ft are shown in Figures 16, 17 and 18. Note that two power 
available lines are shown on each graph. The solid line represents the power 
available predicted by simple theory. The dashed line is a result of the 
ONX/OFFX computer program obtained from Reference (7), and is thought to 
represent a more realistic power available curve. It is clear that the two 
theoretical predictions agree only until approximately M=0.4. With increase in 


speed, the difference between simple theory and ONX/OFFX becomes quite 
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significant. This is important because power available directly relates to excess 
power which in turn is instrumental in defining other performance characteristics 
such as rate of climb and maximum Mach number in level flight. Note also that 


the power required due to drag divergence is not included in this analysis. 
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Figure 16. Power Available and Power Required at Sea Level 
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Figure 17. Power Available and Power Required at 15000 Feet 
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Figure 18. Power Available and Power Required at 35000 Feet 
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D. Climb Performance 

AEW Rate of Climb at sea level and 15000 feet is shown in Figure 19. Rate 
of Climb plots were generated at various altitudes until a service ceiling (rate of 
climb < 100 fom) was found. A plot of the climb rates vs. altitude is presented in 
Figure 20. It was determined the AEW aircraft will have a service ceiling of 
approximately 38260 ft. Although a service ceiling was not specified in the 
Proposed RFP, this ceiling is sufficient to perform the AEW mission. It is 
approximately 1660 feet higher than the service ceiling of the E-2C. Also note 


that the AEW aircraft has an absolute ceiling of 38600 feet. 
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Figure 19. AEW Climb Performance at Sea Level and 15000 Feet 
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Figure 20. Absolute and Service Ceiling Determination 


E. Range and Endurance 

Range and Endurance predictions are shown in Figures 21 and 22 
respectively. Both predictions are made using the Breguet equations obtained 
from Reference (19). The Range and Endurance plots are shown with variation 


in velocity at 35000 ft. 
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Figure 21. 


1200 


1 : : : : : 
: : : Ë : : 
` - . z š 
À : w Е Б ° 
Р : г : : : : 
Е : : : : : 
: • . в ы a 
= — vo meenam е-е * E 
Олин то > ур rere ree + G A Me oe ee ee ва 2 
. • 
: : ; 
: : 


en oo omame a i am qun ——— n 


200 


16 
14 
(2 
lob 
8 

6 

4 

2 

0 


(gx eosuxm pug 


Velocity (fps) 


Figure 22. AEW Endurance at 35000 Feet 
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F. ACCURACY OF PERFORMANCE ANALYSIS 

As with any analysis, it is important to examine the results of the 
performance analysis based on past experience and on historical trends of 
Similar aircraft. In other words, “Are the results of this analysis reasonable?” 

Based on historical trends of aircraft performance, it is clear that the climb 
performance (Figure 19) is far too optimistic. Based on the described design of 
the AEW aircraft, it is very unlikely that it would be capable of climbing at nearly 
12000 fpm at sea level. One possible explanation for this performance is too 
large a T/W ratio. It is unlikely however, that this is a significant part of the 
problem. According to this analysis, even if the AEW aircraft's T/W ratio was half 
the current ratio of 0.46, the aircraft would still climb at sea level at 6000 fpm. 
This is clearly unreasonable. Two other possible explanations of the optimistic 
climb performance are immediately apparent. First, the predicted CDo of may 
be far too optimistic. The CDo analysis does not account for interference drag. 
As a result, the actual CDo is usually higher than the predicted value. This 
difference might be significant on the AEW aircraft which probably has 
Substantial interference drag. It should be noted that the CDo of the E-2C is 
0.0375 which is far higher than the predicted AEW CDo of 0.0205. Second, the 
actual lifting efficiency may be lower than the preliminary estimation. A more 
accurate analysis of the aircraft's aerodynamic characteristics wil! be possible 
only after Computational Fluid Dynamics (CFD) analyses, or wind tunnel tests 


are performed. 
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The results ot the Range and Endurance analyses (Figure 21 and 22) are 
also unreasonably optimistic. Because both the fuel capacity (14000 Ibs.) and 
the TSFC (0.33) are reasonable, it is likely that the aforementioned 


explanations would account for the unrealistic range and endurance results. 
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VI. STABILITY AND CONTROL 


In order to understand what the handling qualities of the AEW aircraft might 
be, a stability and control analysis of the aircraft is necessary. The purpose of 
this chapter is to provide a conceptual analysis of the stability and control 
characteristics of the aircraft. It is important to note that this analysis is a very 
rough approximation. Some of the parameters are the result of design 
approximations presented in previous chapters. Other parameters are 
impossible to predict accurately without the use of wind tunnel testing. In these 
cases, the value of the parameter was selected based on similar existing aircraft 
and past experience. 

The analysis was performed at three mission-relatable flight conditions. 
The flight conditions are: 1) M= 0.2 at sea level, 2) M = 0.48 at 35000 feet and 3) 


M = 0.76 at 35000 feet. 


A. STABILITY AND CONTROL DERIVATIVES 

The stability and control derivative analysis was performed in accordance 
with References (8), (18) and (20). A stability and control computer program 
was written in MATLAB and is included as Appendix |. The analysis assumes 
no aeroelastic effects of the aircraft. All derivatives have the units of rad:!. 


Finally, any effects of thrust have been neglected in this analysis. The stability 
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and control derivatives for the AEW aircraft are shown in Table 9, along with an 


E-2C comparison at M=0.4 and 30000 feet. 


B. DYNAMIC ANALYSIS 

The dynamic analysis was performed in accordance with Reference (20). 
A dynamic modes computer program was written in MATLAB and is included as 
Appendix J. The analysis assumes small perturbation, linear theory. Results for 
the Short Period and Phugoid (or Long Period) modes are approximated to 
second-order systems. Any effects of thrust have been neglected in this 
analysis. The dynamic modes for the AEW aircraft are shown in Table 10. 

The short period natural frequency (Wn) and damping ratio (Z) are 


approximated in Reference (20) as: 
М/п-Ү((2-.“Ма)/и:)-М-)) (1) 
Z--(Mq-«M(. dot)«Z. /ug) (2" Wn) (2) 


A representative example of the dynamic modes is graphically presented in 
Figure 23. The figure shows the short period mode at the three flight 
conditions. All three primary modes have similar characteristics. They are all 


relatively lightly damped with very long periods and small amplitudes. 
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TABLE 9. AEW STABILITY AND CONTROL DERIVATIVES 


DERIVATIVE M=0.2 at M=0.48 at M=0.76 at E-2C 
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TABLE 10. AEW DYNAMIC CHARACTERISTICS 
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Figure 23. Short Period Response 


C. ACCURACY OF STABILITY AND CONTROL ANALYSIS 

One of the advantages of the dynamic analysis is that the final results (i.e., 
damping frequency and period) are directly relatable, and easily 
understandable, handling characteristics. The accuracy of these characteristics 
can be qualitatively evaluated based on historical trends and past experience. 
The accuracy of the dynamic characteristics are directly related to the accuracy 
of the stability and control derivatives, because the derivatives are used in the 
dynamic analysis. 

The results of the dynamic analysis are clearly unreasonable. The most 
obvious discrepancy is in the periods of the three primary dynamic modes (short 
period, long period, and dutch roll). Short period and dutch roll periods for an 


aircraft of this kind typically range from 2 to 8 seconds. Obviously, values 
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ranging between 40 and 206 seconds are unreasonably large. The long period 
values between 1595 and 8770 seconds are also unreasonably large. Long 
period values for an aircraft of this kind are typically about 120 seconds. Also 
note the very lightly damped frequencies of all three primary dynamic modes. It 
is unreasonable that these modes would be so lightly damped, and is 
inconsistent with historical trends. 

Many of the stability and control derivatives appear unreasonable as 


compared with the E-2C. The most unrealistic AEW derivatives include Cm., 
CL(. dot), Cm(. dot), Cmq, and Clp. This would naturally cause unreasonable 


dynamic results. The short period approximation equations are shown on 


page 50. Since Cm. and Cmq are inaccurate, this will result in an unrealistic 
natural frequency. Also, since Cm(. dot) and natural frequency are inaccurate, 


this causes an unrealistic damping ratio. Poor initial assumptions are the most 
likely cause of the unrealistic derivatives. Some inputs were impossible to 
accurately predict within the scope of this research. Such inputs include the 
downwash gradient at the horizontal tail, Cmo, and the moments of inertia. One 
primary conclusion can be drawn from this analysis. Although the method for 
attaining stability and control derivatives in Reference (18) is extremely 
detailed, truly accurate stability and control derivatives can only be acquired 
from wind tunnel tests on a scaled model. Because most of the unrealistic 
derivatives are longitudinally related, any follow-on research should include a 


thorough re-examination of the longitudinal analysis. 
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VII. CONCLUSIONS 


A. ACCURACY 

Because this thesis presents the results of a conceptual design, the 
aircraft's characteristics are by their very nature, a first iteration only. Future 
studies of the AEW aircraft must necessarily include wind tunnel tests of a 
scaled model. Reasonably accurate values of many of the aircraft's parameters 
can only be obtained through wind tunnel tests. 

One of the genuine benefits of this research was the many computer 
programs that were generated. As the design process for this (or any other) 
aircraft continues, these programs can be used to obtain more accurate results 


through the input of more accurate parameters. 


B. EXISTING  ROTODOME/AVIONICS 

Before the design of this aircraft proceeds beyond the preliminary design 
stage, consideration must be given to the use of new airborne detection 
technologies. Based on historical trends, it is likely that the integration of the 
E-2C's detection system into a new airframe will be difficult. The result would 
be an increase in both developmental and life cycle costs. Although new 
detection technologies such as a phased-array radar may be costly to develop, 


the benefits and the life cycle costs must be investigated. 
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C. SUPERCRITICAL AIRFOIL 

Use of supercritical airfoils on aircraft is a relatively new technology that 
should be explored further. The airfoil appears to be ideally suited for aircraft 
that must operate in the transonic regime, and display aggressive endurance 


characteristics. 


D. POSSIBLE PROBLEM AREAS 
1. Escape System 

Within the scope of this design effort, no satisfactory ejection system 
could be determined. The obvious hinderance to a viable ejection system is 
use of the existing rotodome antenna. Difficulties in developing a viable 
ejection system will most likely occur, regardless of the system, as long as a 
conventional rotodome antenna is used. A conventional early warning phased- 
array radar system for example, would be approximately the same size as the 
current antenna. The difficulties in ejection therefore, would be similar. Ejection 
of the aircrew would be much more successful with an antenna that is not in the 
form of a rotodome but within the wings and body of the aircraft. This would 
necessitate the use of a phased-array radar system, and therefore, would be 
costlier to develop. Before a formal AEW HFP is developed, a clear decision 
will have to be made on the aircrew escape system issue, and the resulting 


impact on the radar system. 
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2. Divergent Drag Mach Number (Maa) 

Although the wing Maa of 0.81 is high enougn to operate in the required 
regime, future studies should include an analysis of the drag penalties of other 
aircraft parts in this transonic range. Emphasis snould be placed on the 
fuselage and the rotodome antenna. The relatively wide fuselage and blunt 
nose may cause significant drag penalties at the target high-speed dash Mach 
number of 0.78. With a thickness ratio of 0.3, the rotodome antenna is also 
likely to have a Mgg far below the required operating range. It may, of course, 
require transonic wind tunnel tests to verify how significant these drag penalties 
are. 

3. Horizontal Tail Effectiveness 

It can be seen from Figure 6, that the horizontal tail is directly behind 
the wing and rotodome support pylon. The aerodynamic disturbance created 
by the wing and pylon could result in the loss of horizontal tail effectiveness 
under some flight conditions. This can only be verified however with wind 
tunnel tests of a scaled model, or by a CFD analysis. 

4. Wingfold System 

Another area of difficulty could be in the wingfold system. Because a 
double-wingfold system is new technology, developmental costs may be high. 
The double-wingfold will be an engineering challenge to both the structures 
and the flight control design teams. It should be pointed out that if an aircraft 


design employs a phased-array radar system with a non-conventional antenna 


SN 


such as the one previously mentioned. the need for a double-wingfold system 


might be eliminated. 


E. RECOMMENDATIONS 

Within the scope of this research, the design of an AEW aircraft using the 
existing rotodome and avionics should be abandoned. Use of the rotodome will 
negatively affect the aircraft's normal and emergency operations. Considering 
all factors involved, it is unlikely there will be substantial savings using the 
existing rotodome and avionics. 

Future aircraft designs should include integration of a phased-array radar 
system. This system offers the flexibility needed for an aircraft required to 
possess ejection and wingfold systems. Reference (21) provides an example of 
such a design. The aircraft, called the Boeing EX, is shown in Figure 24. A 
comparative analysis of the Boeing EX and the AEW aircraft is provided in 
Table 11. It is clear from the Figure 24, that the phased-array radar system 
allows for more flexibility in the design process, and eliminates the 


aforementioned ejection and wingfold problems. 
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Figure 24. Boeing EX [Ref. 21] 


TABLE 11. AIRCRAFT COMPARISON 


x CHARACTERISTIC | BOENGEX |  AEWAIRCRAFT | 
Takeoff Weight — 
x TW 
[Antenna 
[Ejection Capability — 











Takeoff Weight 
TW 
Antenna Mounted in Wings 
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In conclusion, it must again be emphasizeq that this analysis was the first 
iteration on a conceptual design only. Therefore, the scope of the research was 
limited. A more complete analysis is only possible after an entire design team 


is assembled. 
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APPENDIX A 


AEW AIRCRAFT DESICH 
NAVAL POSTGRADUATE SCHOOL ' 


PROTECE OBJECTIVES 


The object of this design study is to perform the necessary trade 
studies required to define the most cost effective, low risk 
airframe configuration capable of meeting future airborne early 
warning (AEW) requirements in the 21st century. The mission is a 
deck-launched high speed dash, low speed loiter at 20,000 to 35,000 
feet altitude and return. The goal is to select the greatest high 
speed dash Mach number consistent with the maximum range and loitet 
requirements that will provide a carrier suitable aircraft. The 
aircraft will have ejection capability provisions for all members 
of the four to six member aircrew. A fanjyet (no turboprops) power- 
plant will provide aircraft propulsion. The EX configuration must 
exhibit low initial purchase cost and low life-cycle cost. 


6 1 


MISSION DEFINITION 


DECK LAUNCHED SURVEILLANCE: The total mission cycle time (quadruple 


cycle) 


is desired to be at least 7 hours 30 minutes (with one re- 


fueling) plus reserves with a minimum acceptable cycle time (triple 
cycle) of 5 hours 45 minutes (no refueling) plus reserves. 


l. 


10. 


11. 


For taxi, warmup, takeoff and acceleration to M=0.3; fuel 
allowance at sea level static thrust is equal to 5 
minutes at intermediate thrust (no afterburner). 


Acceleration: Maximum power acceleration from H-0.3 to 
best rate of climb speed at sea level. 


Climb: Best rate of climb to optimum cruise altitude for 
design cruise Mach number. 


Cruise: Cruise-out (High speed dash at M=0.7-0.85) at 
design Mach number at optimum cruise altitude. 


Turn: 3g sustained desired; 2g sustained minimum at the 
weight corresponding to the end of crujse-out. 


Loiter: Conduct surveillance at maximum endurance flight 
condition for minimum of 4 hours 30 minutes (200 mnm 
station, no refueling). 


Descent: Descend to best return cruise altitude (no time, 
distance or fuel used allowances). 


Cruise-back at optimum altitude and best cruise Mach 
number. 


Descent: Descend to sea level (no time, distance or fuel 
used allowances). 


Land. 


Reserves: Fuel allowance equal to 20 minutes loiter at 
sea level at speed for maximum endurance plus 5% of 


initial total fuel. 
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WEIGHT: 


CREW: 


AVIONICS: 


Drie DEPRNSE ? 


LOAD FACTOR: 


CARRIER 
SUITABILITY: 


DESIGN CRITERIA 
The maximum takeoff gross weight will be 60,000 lb,. 


The aircraft will have an aircrew of from four to 
six members, including a single pilot. A weight 
allowance of 230 lb, is required for crew members 
and his/her equipment. 


Design an optimal configuration of flat pànel dis- 
plays for tactical cockpit operation. Nominal dis- 
play sizes for consideration are 6x8, 8x8, 13x13, 
3x5, 6x6 and 4x4. Determine any other feasible 
sizes. Architecture for the operation of the dis- 
plays should not be of concern. Recommend (trade 
study result) the best possible combination of 
displays based on the need for the pilot to control 
the aircraft during takeoff, landing and on-station 
flight; consider also the best display combinations 
based on viewing and interactions with tactical 


displays. 


Data/graphics displayed on a panel of any given 
size should be interchangeable with any other pane! 
of the same size. Consideration must be given to 
supportability (e.g. availability of display sizes 
in other aircraft communities) and to minimizing 
clutter. Recommend screen formats for the transfer 
of as many discrete functions and indicators as 
possible to flat panel displays. Use the existing 
24 foot rotodome. 


Presume that a future missile would be the size of 
a compressed carriage AIM-7 Sparrow and would weigh 
500 lb,. Two missiles are required. A chaff and 
flare launcher is required. Provide two wet wing 


stations. 
3g sustained is desired; 2q sustained minimum at 


the weight corresponding to the end of cruise-out. 


Compatibility with CVN-68 carriers and subsequent 
implies the following criteria: 


1. MK-7 mod 3 arresting gear. 


2; C13-1 catapults. 

3: 130,000 1b, maximum elevator capacity (aircraft 
plus loading plus GFE). 

4. 85x52 foot elevator dimensions. 

DH 57 feet 8 inches minimum station "o" to JRD 


hinge for MK-7 JBD locations. 
6. 18 feet 9 inches minimum from tailpipe to JBD 
hinge. 
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LAUHCII: 


ARREST: 


WAVE-OFF: 


POWER PLANT: 


COCKPIT: 


IN-FLIGIIT 
REFUELING: 


STRUCTURE: 


SELF-DEFENSE 
CAPABILITY: 


GROWTII: 


COST: 


GENERAL: 


du Maximum, unfolded span of 82 feet. 

22 foot maximum landing gear width. 

9. 25 foot maximum hanger deck height except 
under VAST stations in the forward part of the 
hanger where the clearance 15 17 feet 6 
inches. The maximum folded height of the 
aircraft should not exceed 18.5 feet. 


со 


Launch wind-over-deck (WOD) should not exceed zero 
knots operational. Operational is minimim plus 15 
knots. Assume a 5 knot improvement on the С13-1 


catapult. 


Arresting WOD should not exceed zero knots. Assume 
a 5 knot improvement on the MK-7 mod 3 arresting 
gear. Approach speed for WOD calculations is 1.05 
times V approved. 


For multi-engine aircraft, a minimum wave-off rate 
of climb of 500 feet per minute, with one engine 
inoperative, shall be available. 


Fan jets (perhaps, upgraded TF-34 engines). No 
TURBOPROPS. 


Nigh visibility cockpit is required for pattern 
work at ship. 


The aircraft must have an in-flight refueling 
capability. 


The airframe structure must accommodate BIRST. 


The EX aircraft must have a self-defense capability 
[derived from complete (survivability,  vulner- 
ability and susceptibility) studies]. 


The structure must be capable of considerable 
weight growth beyond the initial production 
configuration (at least 4,000 lb,). 


Low purchase cost and low life-cycle cost is highly 
desirable. Assume a total buy of 50 aircraft. 


Attention shall be given to quality, maintain- 


ability, manufacturability and concurrent 
engineering issues. 
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APPENDIX B 


ХТћ!з is a constraint analysis program which Is designed to plot various flight 
Xcondit!ons as a function of thħhrust-to-seight ratio (Tst/Nto) and seing loading 
X(Uto/S).Thie program Incorporates different casss which corresponds to 
Xdifferent flight conXditions. Each case wlll bs seperated with a dashed line. 
Xthis program Ís based on the materlal covered In chapter 2 of NattIngly's (et 
Ха!) alrcraft engins design book. All squattons are from Hatt!ingly unless 
Xspscificaliy stated othsrulss. 


XTs1/Hto wlll hencaforth be knoen as TH. Uto/S will bs known as HS. 

XOperative squatton. 

ЖТНИН$=(В/а)*( (4*5/(8В*Н))*(К1*(п*В*Н/(а*5))^2+К2*(п*В*Н/(а*5))+С0о+Я/(а*5))+1/1*9 
/dtX(h*«U^2/(2*go0))) (ean. 2-11) 

XA parabolic drag polar ls assumed. Therefore K2=0 throughout. 


XCase Í:Constant Alt./Speed Crulss. High Speed Dosh @ H=0.78 & h=30K ft. 
Xdh/dt=dU/dt=0. Constant altitude & no accelerat ton. 

nt=tsXnormal g loading 

Ri=O;%Additlonal drag. Aseumed zsro throughout 

к2=0; ХОгад Curva constant 

Bl=0.905;%Welght Fraction 

К11=0.06; ХОгад Curve constant. Obtalned from Hicolal page Е-7. 
P1*2116*.2360:XPressure at 35K ft. 

1170.78: XHach Huabsr 

CDol*.0345;XDrag coefficient at zero IIft (approximate) 
а!"(1.4/2)%р #1172: ХОупаоа | с Pressure 

RRI=0.3106;XDsnsitu ratio at 30K ft. 
aíi*(0.568*0.25*(1.2-H1)^3)*RRI^0.6;XInstalled full thrott!e thrust lapse for a 
high bypass turbofan (sqn. 2-42) 

Ti *1;Xcounter 

for Ц51=20:5:140; Жића range of wing loading 

HSIH(TI)=HS1; 

THI(TI)=(Bl/al)*(K1I1*B1*HSI/ql+K2+CDol/(BI*HSl/ql));Xthe resulting T/U ratio. 
Xeqn 2.12 

11=11+1; Хсоџп,ег 

end 

US1o*q1/Bl*sart(CDot/K11):;XThe minimum U/S for case 1. 

TH1o*(Bl /al) *(KI1*B:*lS10/q1 *K2*CDol/(B1I*llS10/91)) ; 3The minimum T/H for case | 


XCase le: Maximum Endurance @ 35K ft. 

nls*I;Xnorsal g loading 

Ble=0.6;%Welght Fraction 

Ktle=0.045;%Drag Curve constant.Dbtalned from Hicofal раде Е-7. 

Hls=0.45:%tach Numbar 

gle-(1.4/2)*P1*H15^2; XDynanic Pressure 
aíe*(0.568*0.25*(1.2-t11e)^3) *RR1^0.6;X Installed full throttle thrust lapse for a 
high bypass turbofan (sqn. 2-42) 

Tle sXcounter 


бЭ 


fon HSle=20:5:140;Xthe range of wing loading 

HuStel(T1)9HSte: 
THte(T1)-*(Bte/ale)*(Ktte*Bfe*liSte/qte*K2*CDot/(Bte*llSte/afe)) ;Xthe resulting T/I! 
ratlo. eqn 2.12 

TleT1+?:Xcounter 

end 

WStoe=qle/Bletsqrt(CDol/Kite):%The alnimua H/S for case fe 
THlo=(Ble/ale)*(KiletBlethSloe/qletK2+CDol/(BletNStoe/qle)):%The minimum I/H for 
case fe 


Case 2:Constant Speed Climb. This Is a “snapshot” of the climb oniy. Taken at 
Хап assumed TAS=330 fps, H=0.41, RISK ft. w/ an assumed dh/dt of 4000 fpm. 
XdU/dt=0;: 

n291;Xnorsal g loading 

n2-0;Xfidditi!onal drag. Assumed zero throughout 

P2=0.5646%2116.2:;8Pressure at {5K ft. 

U=433;XUelocltu 

dhdt*6?;XRate of Climb (ft/s) 

H2*0.41:XHfach Humber 

B2*0.975;XHelght Fractlon 

K12*0.05;%Drag Curve constant .Obtalned from Hicolal раде Е-7. 
42"(1,4/2)%Р2%1272; ХОупав! с Pressure 

CDo2-0.0345;XD0rag coefficient at zero IIft 

RR2*0.6295;XDensIty ratlo ot 15K ft. 

a2*(0.568*0.25*(1,2-112)^3) *RR2^0.6:XInstalled full throttle thrust lapse for a 
high bypass turbofan (eqn. 2-42) 

T2*1; Xcounter 

for WS2=20:5:140;%the range of wing loading 

HS211(T2) -H92; 

TH2(T2) *(B2/a2)*(K12*B2*152/q2*K2*CDo2/(B2*152/q2) * 1/U*dhdt ) ;Sthe resulting T/I! 
ratlo. eqn 2.14 

12=12+1; Хсоип(ег 

end 

WS20°q2/B2* sqrt (CD02/K12):X%The mintiaum U/S for case 2 
TH20=(B2/a2)*(K12*B2*NS20/q2+K2+CD02/(B2*NS20/q2)+1/U*dhdt);%The minimum T/H for 
case 2 


¥Case 3:Constant Alt./Speed Turn. Sustained g turn. 

Xdh/dt *dU/dt *0 

п3"2;Хпогта! g loading 

R3*0;Xfidditional drag. Assumed zero throughout 
P3-0.4599*2116.2; XPressure at 20K ft. 

B3-0.85;,XHelght Fractlon 

K1320.045; XDrag Curve constant. Obtalned from Hlcolat page E-?. 
K2*0; XDrag Curve constant 

H3=0.46;XHach Humber 

CDo3».03451XDraq coefficient at zero lift 
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43"(1.4/2)%Р3“1372; ХПупатіс Ргеззиге 

RHR3=0.3332;3Dens!tU rot!ñ!o at 20K ft. 
a3=(0.566+0.25%(1.2-M3)^3)*AAR3^0.6:Xinstalled full throttle thrust lapse for a 
high bypass turbofan (eqn. 2-42) 

13= 1; Хсоипі ег 

for WS3=20:5:140;%the range of sing loading 

HS3HM(T3)=US3; 

TH3(13)= (83/03) *(K19%n3*2*B9*HS3/q3+K2*n3+CD03/(B3*°NS3/q3))sXthe resulting I/! 
ratlo. eqn 2.15 

T3*T3+!:¥counter 

end 

WS30°q3/B3* sqrt (CD03/K13);%The minimum W/S for case 3 
1H3o*(B3/a3)*(K13*n3^2*B3*11$30/q3*K2*n3*CDo3/(B3*11530/23) ) ; The minimum T/H for 
case 3 


XCase 4:Horizontal Acceleration 

Xdh/dt=0O;constont oltltude 

n1*1;Xnormal g loading 

81-0; XRddltional drag. Assumed zero throughout 

Ul-100; XInitial velocity. 

Uf *?76;XFIlnal velocity. 

dt*300; Time for acceleration (In seconds) 

P4*2116.4*0,2360; Pressure ot 35K ft. 

dUdt «(Uf -Ul )/dt:XRicceleration 

B1-0.85; XHelght Fractlon 

K11*,055;XDrag Curve constant. Obtalned from Hicolal page Е-7. 
K2-0; XDrag Curve constant 

H4*.58:XHlach Number. A "snapshot" In the alddie of the run 
CDo1-.0345;XDrag coefficient at zero |111 
g*32.17?;XRcceleration due to gravity (ft/sec) 
q4=(1.4/2)*P4*H4“2:XDUunomic Pressure 

AA4*.3106;%Density га о аб 35K ft. 
a4*(0.568*0,25*(1.2-14)^3)*RR4^0.6;XInstalled full throttle thrust lapse for a 
high bypass turbofan (eqn. 2-42) 

2*1/g*dUdt ; 

T4951: Xcounter 

Гог WS4=20:5:140;X%the range of wing loading 

5 411(74) =1154: 

114(74) = (84/04) *(К14+В4+1154/94+К2+С004/(84* 154/4) +2); {Һе resulting T/U ratio. 
ean., 2.18 

T4-*T4*1:Xcount er 

end 


XCase 5: Takeoff Ground Roll 

Xdh/dt *0; 

$g*3000; XGround rol! takeoff distance 
Аћ5 = ,0023769: Хбва leuel депэ {и 
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Kto*l.2;Xstall-to-takeoff velocity ratio 

Cla-2.5;XHax lift coefficient for takeoff 

B5*f!;XHelght Fraction 

H5-«0:XHach Humber 

RR5*!;XDenslty ratio at sea level 
a5*(0.568*0.25*(1.2-05)^3)*RR5^0.6;XInstalled full! throttle thrust lapse for a 
hlgh bypass turbofan (eqn. 2-42) 

9=32.17; ХЯссе | его! Гоп due to gravity (ft/sec) 

T5*1: Xcount er 

for WS5"20:5:140,%the range of wing loading 

HS5H(TS)=HS5; 

THSR(T5)=((20.9*HS5)/(RRS5S*Cl|m))/(Sg-87*sqrt (HS5/(RR5*C|m)));Xthe resulting T/H 
ratio. This !s froa Hicola! (eqn.6-3)1 

T52T5*1;Xcounter 

end 


XCase 7:Landing Rol! 

Xdhdt *0: 

Clm=3.0;XHax I!ft coefficient for landing 

$|25000; XLand!ng distance 

AR=1;%Density ratlo at sea level 

0=0.2:.1:1,2: 

Н58= (51-400) *АА*СІж/ 116; ХЕгом Нісоіа! (aqn. 6-5).Note It Is Independent of T/U. 
for S*1:11, 

HS8n(S) 1S8; 

end 


ХСазе 9: Halntalnabilitg 

HHFHe30; Xltalntenance тап hours per flight hour 

T9s=l;Xcounter 

for WS9*20:5:140,%the range of wing loading 

HS9H(T9)=HS9; 

TH9(T9)«(HfFH/7.25716)-(0.196568/7.25716)*159; Xthe resulting T/H ratio.Thls ts 
XHewberry's equation for the fighter alrcraft only. 

TU9T(T9)- (HHFH/13.6383)-(0.1555/13.6383) *59; the resulting T/H ratio. This !s 
XHewberrg's equatlon using all25 alrcraft. It was used because It ls probably 
Xmost realistic. 

19=19+1 ; Хсоипіег 

end 


plot (HSIN, THI, UWStelt, Tle, S21, H2, 'x' , 53H, H3, ' * , HS4H, T4, '0' , HSSH, THSR, 7**,158 
n,TH8,'-',HuS9n, TH9T, ' -, ") 
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APPENDIX C 


XIhis is on ejection progron sith expressions from iloerner' s Fluld Dunomic Droq 
book, Chopter 13. 


H=3DD;Xselght of the seat nnd creg nenber 

g=32.2;Xoccelerotlon due to gravity 

11*,2:XHach number 

оћН=|.4;Ховата 

P»2116:13*.8321: Xpressure 

q=(GAN/2)*P*t*2:Xdynamic pressure.assumed constant 

Оа=9; Хдгад area (varles beteeen 4 and 9ft^2) 

6=60;Xappraximate average vertical velocity 

Q=1:Xcaunter 

for Y*0:14, 

VH(0)=V; 

T(Q)*V/w;Xt|ne 1s equal ta velocity diuided by distance 
T2(0)*T(Q)^2;Xt Ine squared 

41(0)-8*(g*q*T2(Q)*(0q/U)) ; the front seat trajectary. eqn. 26, chap 13 
K2(0)9#16+(g*q*l2(Q)*(Dq/N)):Xthe back seat trajectary. eqn. 26, chap 13 
0*0*1:Xcounter 

end 


Рат, "Ка, "у, 


Xthls dross the rotodome antenna 
Ru=(9,7413 10.929 9.7413]; 
Rl=(9.7413 9.7413 9.7413]; 
Ac=(9.7413 6.553 9.7413]; 

#0-[16 28 40]; 

plot (XD,Ru,3D,RI, '-', 4D,Re, '- ^), 
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APPENDIX D 


XThis welght program has two parts. The first Is a subroutine which computes the 
Xeelght of thé propulélon and fuel systeas. These figures are needed for the 
Zealn ргодгав which iterates a takeoff welght. 


XThe below values are Inputs that are required for the equations that have been 
Zobtained from “The Fundamentals of Aircraft Design” by Leland N. Hicolla 
(Chapter 20) 

В | =р] *2.375^2; inlet Area 

Ніе2; XHuaber of Inlists 

Kgeo*-!; XDuct Shape Factor 

P2724; Xlax Static Pressure at Engine Compressor Face-psia 

Kte=l: XTeaperature Correctlon Factor 

Каті; Х0исі Haterlal Factor 

Ld-3; XSubsonic Duct Length 

Fge-2154; XTotal Hing Fuel In Gallons 

Fgf=0; XTotal Fuselage Fuel In Gallons 

((-55; %Fuselage Length 

He-2; XHuaber of Engines 

B-72; XHing Span 

Heng-2000; XHelght of Engine 


XThe equation nuabers from Hicolal are included with the appropriate equations. 
Ht fde7.435*Hi *(Ld*Ri^.5*P2)^. 731; 420-15 
Hssc-11.6*((Fge*Fgf)*10^(-2))^.818,120-16 
Hbsc*?.91*((Fgs*Fgf)*10^(-2))^.854,320-18 
Hi fr*13.64*((Fge*Fgf)*10^(-2))^.392;X20-19 
Hdd-7.38*((Fge*Fgf)*10^(-2))^.458; 120-20 
Htp*28.38*((Fgse*Fgf)*10^(-2))^.442; 320-21 
Hec-88.46*((Lf *B)*He*10^(-2))^.294; 320-23 
1as579.33*(He*lleng*10^(-3))^1.078; 120-26 
Hfs=Hssc+ibsc+HddtHtp+HI fr, 

Поре fd*«life*llec*iliss*(Heng*2), 


XThis program ls deelgned to find the appropriate takeoff ив! ОР (Шо) өһеге the 
Xequation Is a polynoalal wlth fraction exponents.The secant method Is used to 
Xfind the desired root.The operative equation (which Is so designated below) 13 
Xset up so that Xthe prograa eli! find lito (a.k.a. X) when Y ls equal to 
Xzero.The many equations that preceed the operative equation are partlons of the 
Xfinal equation. They are seperate to make the operative equation more 
Xaanageable. 


The below valuee are Inputs that are required for the squatlone that have been 


TO 


Xobtalned from “The Fundamentals of Aircraft Design™ by Leland ft. Hicolla 
(Chapter 20) 

Н=4.5; Ultimate Load Factor 

{0с=0. 12; StMaximua Thickness Ratio 

1 1е=(21%р1/180); XLeading Edge Sweep 

Ct*4; XChord Length at Tip 

Сг=13.75; XChord Length at Root 

I*Ct/Cr; XTaper Ratlo 

8-8.11; Xfiepect Ratlo 

66-639; ХИіпа Пгеа 

Sht=180; ЖНогіг?отіа! Та!1! Ріппіога Area 
Bht*24; XSpan of Horizontal Tall 

tRht*0.86; XThickness of Horizontal Tall at Root 
Cmac=9.?77; %HRC of the Uing 

Lt=25; Tall Homent ñra 

Htllu=0; Xiorizonta! Tall Height to Vertical Tall Ilelght Patio 
Sut*45; XUertical Tall firea 

Пе, 768; ХПах!аџа Hach Nuaber at Sea Level 
Сге22; ХАиддег Area 

Rut*«1.111; XRepect Ratlo of Vertical Tall 
|(>0,5; ХТарег Ratlo of Vertical Tall 
Lut*(30*p1/180); XSweep of the Vertical Tall 
q=600; %Haximum Dynamic Pressure 

Lngth*55; XFuselage Length 

H=0; XHaximum Fuselage Helgth 

Кіпі»!; Хілпіей Constant 

Hpll=2; XHumber of Pllots 

He=2; XHumsber of Enginee 

Htron=10000; XHelght of Rulonics 

Hcr*4; XHusber oí Crew 

Ksea*149,12; XEjJectlon Seat Constant 
Hrad*«3086; XRadome lUelght 

Hfuel-*14000; XTotal Fuel Helght 


XTthe equation nuabere from Hicolal are Included with the appropriate equat lons. 
XThe first loop le used to compute the first two valuee of V after the two 
Xinitial guesses for lito (H) haue been made. Two initial guesses are required 
Xfor the secant sethod. 

Pel; 

for Wto=40000:10000:50000,%40K & SOK are the teo Initial guesses. 

Х(Р)= Но; 
le*19,29*(1*H*!ito/toc*((tan(LIe)-(2*(1-1))/(R*(1*1)))^2*1)2*10^(-6))^.464*(CT* 1) *n 
)^.7*9«^.58, 320-2 

Vh*(to*H)^.813*S5ht ^.584* (Bht /tRht )^.033*(C2ac/Lt )^.28; 120-3a 

Uht *.0034*Yh^.915; 120-3a 
Vus(1*HtHu)^.5*(1t0*H)^.363*5ut^1.089*8^,601*L0^(-.726)* (1*8r/S0t )^.21 7*Rot ^. 337* 
(1*1t)^.363*(cos (Lut ))^(-.484) :120-3b 


ү! 


Hut«2*0.19*70^1.014; 1320-3b 
Hf*11.03*(KInl^t.23)*(a*10^(-2))^.245* (0to*10^(-3))^.98* (Lngth/H)^.61;X20-5 
Ш19= 129. 1% (Н10*10^(-3))^. 66; 9320-7? 

Ићуд=2г3, 77"(Ико%107(-3))71.10;Х20-35 

Hf 1*Hpl1*(15*.032*1t0*10^(-3)):120-39 

He1-7He*(4.80*,006*11t0*10^(-3)); 120-40 

Ha1*.15*(Hto*10^(-3)); 420-42 

Пее" 346. 98+ ( (Нее + Ні гоп) *10^(-3) )^. 509: 1320-44 

Це! =Кева+Нсг ^1, 2: 420-50 

Шох" 16.89+*Нсг^ 1, 494; 20-51 

Hac*201.66*((Htron*200*Hcr)*10^(-3))^.735; 120-65 

Hfc*1.08*(1to)^.?;Sthle equation !s from Roskam Port Ü 

XThe below equation Is the operative equation. 

V(P)«*(-Hto)*Hwstiht Hut +H fH i g+Whydsl fl +Nel+Hml +Heatiat *Hox+HactHrad+Wfuel Htronsll 
pp+Hfc; 

PeP«1: 

end 

XThis concludes the loop that computes the values of Y for the two Initial 
Xquesses. 


ХТНе еесопд loop 16 deelgned to actually find the root.The loop allowe for up to 
416 Iterations. 

for J*3:12, 

4(J)*M(J-1) -V(J-1)*((M(J-1) -4(J-2))/(Y(J- 10-V(J-2) D ; XThls 19s the secant method 
Xforaulal tt computes a value of X (Uto) from the previous two K's and thelr 
respective Y valuee. The reet of thle loop juet computes the new value of У 
Хігов the newly corpuXted H. Ноге Information on the eecant method can be found 
Xin any numerical rethods book. 

Hto*4(J); 
He*19.29*(1*H*lito/toc*((tan(Lle)-(2*(1-1))/(R0*(1*1)))^2*1)*10^(-6))^.464* ((1*1)*R 
)^.7*5в^.58:Х20-2 

Vh*(1ito*H)^.813*S5ht ^.584* (Bht /tBht )^.033*(Cnac/Lt )^.28; X20-3a 
Hht*.0034*Yh^,915; 420-3a 
Vus(1*HtHu)^.5*(1to*H)^.363*5ut^1.089*^.601*Lt^(-., 726)* (1*S8r/80t)^.217?*00t^.33?* 
(1*1t)^.363*(coe(Lut ))^(-.484):1320-3b 

lut -2*0.19*V0^1,014:X20-3b 
Uf*11.03*(KInl^t.23)*(q*10^(-2))^.245*(1t0*10^(-3))^.98* (Lngt h/H)^.61; 420-5 
11g-129.1*(1t0*10^(-3))^.66; 320-7 

Hhyd*23.77*(0to*10^(-3))^1.10; 420-35 

ШіеНрі1%(15»,032%!%0%107(-3)):Х%20-39 

He17He*(4.80*.006*1t0*10^(-3)); 120-40 

На! =, 15+ (010+10^(-3) ); 420-42 

Пе=~ 346.96 + ( (НГ + гоп) *10^(-3) )^. 509; 20-44 

Ше =Кева*Нсг^ 1.2, 20-50 

Hox*«16,89*Hcr^1.494:X20-51 

Hac*201.66*((Itron*200*Hcr)*10^(-3))^.735, 320-65 

Hfc«1,08*(Hto)^.7)Sthis eauat!on 18 from Roekam PartU 
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XThe belos equation is the operative equation whos root we are seeking. 
Y(J)*(-Hto)*Hs*tiht Hut «uf *tgeühud*HfT*Hel «Hal *Hes*tist*tlox*Hac*trad*Hfuel*litron*il 
pp+Hfc; 

end 

dlsp(Hto), 

ХИјо= 5.14906+04 Ibs 
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AEW1.XLS 
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MOMENT ARM REFERENCED FROM `S` FEET 


IN FRONT OF THE NOSE, “75 FEET BELOW THE FUS 










WING (OUT) 0 2250 34| 765001 12| 27000 23 
WING (WED — | 3%) 30) 107400) 12] 42960 75 
HORIZONTAL TAIL 445 55.5| 24697.5| 15) 6675 4 56 
NACELLES ` | 969 25.5| 24709.5| 10| 9690 9.75 
FUSELAGE 2757| 29] 79953 9| 24813 O 









— —— — y —— — —- 








VERT TAIL 






FUEL 






2| 168000 





WING 













~" 
Qni 





BLADDER (M) 








оз! 


i 


ООМР5 АМО ОВАК(М) | __ 30) 3) 990| 12| 360 












m 
сл! 





CELL BACKING (M) 






TRANSFERPUMPS(M) | 500 9) 3300) 12] 1320] 75 
INFLIGHT REFUELING | 45] 15 ___675| 10) 450 35 
ENGINES | _ 4000) 25.5| 102000) 10| 40000 975| 
ENGINE CONTROLS | 116] 20| 2320 10| 1160 45 





STARTING SYSTEMS __ 












HYOs — 











LANDING GEAR (NOSE)| 236 -13| 3068] 22] 51921 0 = 
LANDING GEAR (MAIN) | 1473] 39| 57447] 221 32406 6 
HYO SYSTEM |] 1762] 30] 52860 — B| 14096 0 
FLIGHT CONTROL SYS. |  2043| 30] 61290 8| 16344 0 
















Fur INST р 3 10] __ 330] 10] 330] 0 
ENG INST — 7 10 io] 100] 10] 100 0 
AIR CONO | 1159 3i|  35929|  12| 13908 0 
OXY SYSTEM _ | 134 15] 2010 7 938 0 














ELECT SYSTEM | 1165] 35] 40775 .9| 10485 0 
Misc INST) . | 7| 1o 70 ii 7 20 
APU = — Ор 50 25| 1250 ЕЕ 0 
АМОМС5 _______| 10001 41| 410000] 10 0 
RADOME — .— | 3000 33| 99000) 19 0 
CHAFF/FLARE LAUNCH | 300 3 9900 8 0 





82 





AEW1.XLS 


SEATS о м и 8 — 179.9021 
нк | — 


ООО Т М: 
BEEN п — | C | | 
БЕН ООо —| — | | | 
ЕЕЕ: 77513902 EE — 560703] ——— | — — | — 
XCG FROM "5" FEET FORWARD OF NOSE т 0 
ШЕШЕН ЕЕ MEMO. X. EN 
ЭЭ [32.44276] —-— m 
ИШЕ. ШУ = 
ZCG FROM "5 FT BELOW FUSELAGE Ec Ea _ 
EN | 1 | 


| | 





NE^ 

| Б. 
І22-| 14769321 ішіН22 Jo | || 
ББ вы 05ш-/? Ma ` 

Ixz= ЖУ 
12у=| Ol slugs/ft^2 
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1.18786 
1.18786 


_ 16.7272 
0.828301 


(Y-Ycg)^2(2i-2cg)^2| — 


205627.5 


16696.75 


Í — — — = а 


92918.19 


| 10088 97 


244637.8 


- —  — Á У 


47107.51 


42169.57| 32110.6 


39| 33723.89| 177291 


1.18786| 29465.62| 3595.757] 31842.63 


_ 804130 


18786| 1115.636 


1.18786 


25| 0.828301 


m w s w 


0.828301 
| 0.826301 


75.86602 
15.86802 


8.468742 
0| 8.468742 


0.828301 
0.828301 


8260.727 


764.0107| 6765.784 0 


6318.165] 771 


588.5235 


 —— — — = — 


383563.2 
2445.083 


13681.46] 14195.44 


17904.38 | 1068422 
164778.7 


14921.92 


17301.64| 29194.79 


27.33393 
ПА 


1316.729 
2048.721 


4250.527 


0.056561 


2| 423.4371 


3689.968 


42678.09| 40629.37 


12048.81] 7798.287 
3516.379] 3516.323 
3170.888 


сісісіс 
ЗЕЕ: 


| e 


i 


| 


81811 12181813 


278,153 35 


__14.252.7 
-512.68 


0.828301) 8283.008| 745689.8 737 
5.19631] 196338.9] 197373.5)_ 1039.551 
8.468742| 2540.623| 2644.078| 103.4551 


oioio|sioiel 
2|=|= 
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—— 
141583! 0.001 0.00 


ПИШ е Же 
1564.88| 143147.9 
ku АЕ. - да. 


Е-% 

КЕНЕН о ЖЕК wm | !' | 
И = 
"s. —— | 

3217604| 2386534| 4751882| 0.00] 0.00 


100006.3| 74175.85| 147693.2 ШЕҢ 
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APPENDIX F 


Coordinates of 12 Percent Thick Supereritieadl Virfoil SC(O) 0712 
Designed (от 07 110. Coetlicient 
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APPENDIX G 


X2ero llft drog coefficent of entire aircraft. This progroa alll compute 
Xlsoloted ports of the olrcroft & then sum them. This 13 from ORTCON. 


XPart 1: lsoloted Шпа 

Cr*13.75;XRoot Chord (ft) 

Сі-4 ХТІір Chord (ft) ' 
toc*.12;XThlcknese Ratlo 

Lle=2!*p1/180;%Leoding Edge Sweep (rads) 

8=72: ХИ па Span (ft) 

HU=1.573*10“(-4);XUlscositu (ft“2/s) 

Vinf=820;XFreestream Velocity (ft/s) 

|=Ct/Cr;XToper Rotlo 

B2=B/2; Half Hing Spon (ft) 

TLle-ton(Lle);XTangent of Leoding Edge Sweep (rods) 

Ctp*TLIe*B2; 

Crp*Ct*Ctp-Cr; 

6fp*«2*((B2*(Cr*Crp))-(.5*B2*Ctp)-(.5*B2*Crp)); Xing firea (ft^2) 
Cb*(2/3)*Cr*((1*1*1^2)/(1*1)):XC bor - leon Rerodunamlc Chord 
Re-UInf*Cb/HU;XReynolds Huaber 
Cbf*0.455*(10gl0(R&))^(-2.58); XRveroge Turbulent Skin Friction Coefficient 
Cdos-2*Cbf*(1*(2*toc)*(100*toc^4)),XCdo of the ling. ean. 4.1.5.10 


XPort 2: Isolated Rotodome (not Including Pylon) 

Сгр=24; ХЛободове Root Chord (ft) 

Ctr*0:XRotodoae Тір Chord (ft) 

tocr*.135; XRotodoae Thickness Rot lo 

Ir*«Ctr/Crr;XRotodoae Toper Rotlo 

Cbr*(2/3)*Crr*((f*lr*1Ir^2)/(1*1r)); XC bor - Rotodome Neon Rerodynamic Chord 
Rer*Uinf*Cbr/HU; XReunolds Humber 

Chfr=0.455*( log!O(Rer) )*(-2.58);Rotodome Average Turbulent Skin Friction 
{Coefficient 

Cdor*2*Cbfr*(1*(2*tocr)*(100*tocr^4)):XCdo of Rotodome prior to multiplication 
Хо! Rotodose-lling Ягео Ао! Іо. едп. 4.1.5.10 

Sr=p{*12^2:%Rotodome Area (ft^2) 

Cdorp*Cdor*Sr/Sfp,XCdo priae of Rotodome 


XPort 3: Rotodome Pylon (Support) 

XThe Pylon has been approxiaoted oe a wing with the following dimensions. 
Crs=!3;X%Rotodome Pylon Root Chord (ft) 

Cts=8;XRotodome Pylon Tip Chord (ft) 

tocs=.3;%Aotodome Pylon Thickness Rotlo 

le*Cts/Crs;XRotodose Pylon Taper Rat lo 

СЬз=(2/3)*Сгз*((1+1я+1#^2)/(1+1е)),Ж С bor-Rotodome Pylon Меоп Aerodynamic Chord 
Res=U|nf*Cbs/HU;XRéunolds Humber 

Cbhfs=0.455*( foglO(Res))*(-2.58):¥Rotodome Pylon Average Turbulent Skin Friction 
Coefficient 

Cdos*2*Cbfs*(1*(2*tocs)*(100*toca^4))1XCdo of Rotodome Pulon prion to 
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Xmultipilcation of Pulon-Hing Rrea Ratlo. eqn. 1.1.5.1a 
$97((13*:0)/2)*6.4; XRotodone Pylon Area (ft^2) 
Cdosp*Cdos*Ss/Sfp,XCdo prime of Rotodome Pylon 


XPart 4: leolated Fuselage (Body) 

XThis program assumes a oglve shaped body. 

Овох=8; ХНах Dlameter of Fuseiage 

(бе55; ХРизе! аде Length 

FA=Lb/Daax:%Filneness Ratlo 

Db*1.0;XBase Diameter 

Rebs*UInf*Lb/HU; XReynolds Husaber 

Cbfíb*D.455*(logl0(Reb))^(-2.58); XFuselage fiverage Turbulent Skln Friction 
XCoefflclent 

SwoSb916.65;3From USAF S&C DatCom Figure 2.3.3 

Sb*pl*4^2;XFrontal Area of Fuselage 
Cdof=t.D2*Cbf*(1+(1.5/(Lb/Dmax)^1.5)+(7/(Lb/Dmax)^3))*SwoSb; %Cdo-Fuselaga Skin 
XFrictlon. First part of eqn. 4.2.3.la 
Cdobb*(D.029*(0b/08ax)^3)/(sqrt(Cdof)); XBase Pressure Cdo. eqn. 4.2.3.1b 
Cdob-Cdof*Cdobb;XCdo of Fuselage prior to multiplication of Fuselage-Wing Area 
Хћа| | о. едп. 4.2.3.!а 

Cdobp=Cdob*Sb/Sfp,%Cdo prime of Fuselage 


Part 5: Isolated Horlzontal Tall 

Crh=9;XHorl zontal Tal! Root Chord (ft) 

Cth=6;%Horlzontal Tall Tlp Chord (ft) 

Cthp=3; 

toch=.12;:fHorlzonta! Tall Thickness fatlo 

Bh2*12:XHorlzontal Tall Half Span 

lheCth/Crh;XHorlzontal Tall Taper Ratlo 

Cbh*(2/3)*Crh*((1*Ih*Ih^2)/(I*1h)); XC ber-Horl zontal Tall Mean Aerodynamic Chord 
fieh-UInf *Cbh/HU; XReunolds Humber 

Cbfh*0.455*(log10(Reh))^(-2.58);XHorlzontal Tall Ruerage Turbulent Skin Friction 
XCoefficlent 

Cdoh-2*Cbfh*(1*(2*toch)*(100*toch^4));$Cdo of Horizontal Tall prior to 
Xmultipiication of Horizontal Tall-Wing Area Ratio. eqn. 4.3.3.1a 
Saph*2*(Crh*Bh2-,5*Bh2*Cthp);XHorlzontal Tall Rrea (ft^2) 
CdohpeCdoh*Saph/Sfp,XCdo prime of Honrlzontal Tall 


XPart 6: Isolated Uert!ical Tal! 

Cru*6:XUertical Tal! Root Chord (ft) 

Ctuv*3; XUertlcal Tall Т!р Chord (ft) 

Cthp*3; 

tocu*,12;XUertical Tall Thickness Ratlo 

IvsCtu/Cru;XUertical Tall Taper Ratlo 

Cbu*(2/3)*Cru*((1*1u*1u^2)/(1*10)) :4C bar-Uertlcal Tal! Mean Aerodyunamle Chord 
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Reu*Uinf*Cbu/HU; XReunolds Humber 

Cbfu*0.455*(logl0(Rev))^(-2.58); XUertlcal Tall Ruerage Turbulent 5Кіп Frlict lon 
XCoefficient 

Cdous2*Cbfu*(1*(2*tocu)*(100*tocu^4)) :XCdo of Uertical Tal! prior to 
Xmultiplication of Vertical Tall-Wing Area Ratlo. eqn. 4.4.3.1a 
Sapu*90:XUertical Tall Area (ft^2) 

Cdoup*Cdou*Sapu/Sfp,XCdo prime of Vertical Tall 


XTotal ; 
Cdo=Cdow+Cdorp+Cdosp+Cdobp+Cdohp+Cdoup,XTotol Alreraft Cdo. ean.4.5.3.1b 
Cdoa=Cdow+.008+Cdobp+Cdohp+Cdoup, Xfotal Rircraft Cdo using actual rotodome drag 
Information. 


XCdo =0.0177 
XCdoa*0.0205 
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APPENDIX H 


Xthis program is designed to calculate the Coefficient of Orag, LIft-to-Drag 
XRatio, Thrust Required, Power Required, Power Avallable, Excess Poser, Rate 
Zof Climb, Endurance and Range. Thé equations are found In any !Intrductory 
Xalrcraft book. This analXysis was performed using Anderson's "Introduction to 
XFilght, Chapter 6. 


Cdo=0.0205;%Alrcraft Coefficient of Orag 
RR-8.11;XRspect Ratlo 
e*0.8; XEfficlencu 
1*53000; XRircraft Welght 
Ufuet-14000; XFuel Helght 
1e*53000-14000; XEnptu Це! ана 
n0-.0023769*1; X0ensitu (s1/ft^3) 
S1G-R0/.0023769;: XDenslItu Ratlo 
Thr*25400* (S16) ; XThrust 
SFC*0.33/3600; XSpecific Fuel Consumpt lon 
5639; ХИ па Агеа (1472) 
K*1/(pl*RR*e); 
T=1:Xcounter 
for R=.05:.05:3,XThls ts the range of Cl chosen. 
CI(T)-R;XCoefflclent of Lift flatrix 
Clsq(T) R^2; XCI1 squared 
Cd(T)=Cdo+K*A*2;XComputed Cd Matrix. eqn. 6.1c 
LoO(T)-CI (T)/Cd(T) ; XLI ft-to-Drag Ratlo (max L/0=16) 
TR(T)*U/LoO(T);:XThrust Required for Level, Unaccelerated Flight. ean. 6.15 
U(T)*sqrt(2*U/(RO*S*CI(T)) ) ; Ue loclty calculated from Ci. eqn. 6.16 
PTACT)=.5*AO*U(T)*2*S*Cdo;8Porasitic Thrust Required for Level, Unaccelerated 
XFlight. eqn. 6.17 (tet part) 
ITR(T)*,5*RO*U(T)^2*S*k*R^2; XInduced Thrust Required for Level, Unaccelerated 
XFlight. eqn. 6.17 (2nd part) 
PR(T)-TR(T)*U(T) ; Power Required for Level, Unaccelerated Flight. eqn. 6.23 
PRp(T)*sqrt(2*1^3*Cg(T)^2/(R0*S*CI(T)^3)) ; SPower Requlred for Level, 
XUnaccelerated Flilght (double check). eqn. 6.26 
PPR(T)-PTR(T)*U(T) ;SParasitic Power Required for Level, Unaccelerated Flight 
IPR(T)- 1TR(T) *U(T) : Si nduced Poser Required for Level, Unaccelerated Filght 
PRp(T)-Thr*U(T):XPower Rualloble (the slope of thls line Is the thrust) 
EUR(T) *(1/SFC)*LoD(T)*1og(H/lle) ; S€£ndurance. eqn.(65.63) 
RHG(T)*2*sqrt (2/(R0*S))*(1/SFC)* (sqrt (CE CT))/Cd( T)) *(sart (H) -sart (He) ) ; Range- 
Xeqn. (6.68) 
Gang(T)*atan( 1/LoD(T))*(180/pl);XGl Ide angle (In degrees). eqn. 6.47 
XGrng(T)=H*LoO(T);XGitde Range. figure 6.30 
T*T*1:Xcounter 
end 
4*1: Xcounter 
for Ufi*0:35.7:999,6, X0 to 1000 fpe 
UR (W) UR; XUelocitu Matrix 
PA(X)=Thr*UA: Power Avallable Hatrix (Thr Is the slope of thls line) 
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ХеХх+ |; Хсоип(ег 
end 
Рб=РАр-РА; ХЕхсезз Ровег Hotrix. eqn. 6.42 
RoC=PS/U:fAate of Climb. eqn. 6.43 
Thet=asin(Aol./U).*(160/pl):%climb angle. eqn. 6.41 
Xdisp(LoD), , 
Xdisp(PS), 
Xd!sp(RoC.*60), 
Yplot(Cd,Cl), 
Xplot(Cd,Clsq), 


XYplot(U, TR), 
Xplot (U, TR,U,PTR, '--,U, ITR, '--"), 
Xplot (U,PR,U,PPR,'--',U,IPR, '--',U,PRp, x ,URH,PR, ' -), 


Xplot (U,EOR./3600), 
Xplot(U,RN6./6000), 
Xplot(U,RoC*60), 


Ythis is a result of actual thrust/poser obtained from OHX/OFFX 

РЯ=1=[0347933 11130570 13370120 13693171 14048422 13970359 13052273]; Хасіџаі РП 
Hatrix at sea level 

PAIS =1.0e+07*(0.5347 0.70064 0.6346 1.13623 1.2283);%Power Avaliable at 15K 
PAIS =1.0e+06*(2.2604 3.0139 3.6222 5.5050 6.2335] :%Powser Avaliable at 35K 
Hsl=[,3 ,4 ,5 .6 .7 .8 .9): 

H=U,/(1116); 

Пат" ЈАН. / (1116); 

HIS=[.3 .4 5.009): 

135=[,3 ,4 ,5 ,8 .9): 

Хр!о (И,РА, 'о' ,Н,РАр, ' -' „Лам,РА, '"-' ,И,РА, ' -' ,135,РА35, ' --'), 

PSRI*1.0e*07?7*[0 .2195122 .6585366 .8780488 1.0341463 1.03 .9993 .9405 .8893 
‚9443 .9046 .7692 .7374 .7087? .6825 .6586 .6365 .6161 .5972 .5796 . 5631 .5477 
.5332 .5194 .5065 .4942 .4825 .4714 .4609 .4508 .441!1 .4318 .4230 .4144 .41062 
.3983 zeros(1,25)]; 

PSR2-*1.0e*07*[zeros(1,36) .3907? .3834 .3763 .3694 .3628 .3563 .3501! .344! .3382 
‚3325 .3269 .3215 .3163 .311 .3062 .3013 .2966 .2920 .2874 .2030 . 2787 .2745 
‚2703 .2663 .2623};: 

PSR=PSRI+PSA2; factual PS (excess power) matrix at Sea Level 

НА! + (,04 ,8 .7 .6 .5 .45 .4198 .3086 .3635 .3427 .3252 .3100 .2968 .2852 .2748 
.2655 .2571 .2494 .2424 .2359 .2299 .2244 .2192 .2144 .2099 .2056 .2017 .1979 
.1943 .1909 .1877? .1847 .1818 .1790 .1763 .1738 .1714 .1690 .1!668 .1647 .1626 
zeros(1,20)]; 

ИА2=[260го9(1,41) , 1606 . 1507 .1568 .1550 .1533 .1516 .1500 .1404 .1469 .1454 
.1440 .1426 .1412 .1399 .1386 .1374 .1362 .1350 .1339 .1327]; 

ИА=ИЯ1+ИЯ2, 

RoCA=(PSA./U)*60; factual Aol Matrix 

Xolot (HA, AoCA), 

PSAISt=1.00+06*(0 1.852 4.259 5,556 6.204 6.296 5.926 5.6713 5.4431 5.2319 
5.0362 4.8543 4.6846 4.5260 4.3771 4.2371 4.1051 3.9804 3.6621 3.7499 3.6431 


$1 


3.5413 3.4441 3.3511 3.2620 3.1766 3.0944 3.0154 2.9394 2.8660 2.7951 
zeros(1,26)]; 

РбЯ152=1 ,0е+06* [хегоз(1,31) 2.7267 2.6605 2.5963 2.5342 2.4739 2.4154 2.3585 
2.3032 2.2494 2.1970 2.1460 2.0962 2,0477 2.0003 1.9541 1.9089 1.8647 1.8215 
1.7792 1,7378 1.6973 1.6575 1.6186 1.5804 1.5430 1.5062]: 
PSRIS«PSRISI*PSRIS2; Xactual PS (excess poser) matrix at 15K 

HAISt-{.957 .9 .8.7 .6 .5 .45 0.4023 0.38652 0.3701 0.3566 0.3445 0.3336 0.3236 
0.3145 0.3061 0.2984 0.2912 0.2845 0.2782 0.2724 0.2669 0.2617? 0.2568 0.2522 
0.2478 0.2436 0.2396 0.2359 0.2323 0.2288 0.2255 0.2224 0.2194 0.2165 
zeros(1,22)]; | 
В152={гегоз(1,35) 0.2137 0.2110 0.2084 0.2059 0.2035 0.2012 0.1989 0.1967 
0.1946 0.1926 0.1906 0.1887 0.1868 0.1850 0.1833 0.1816 0.1799 0.1783 0.1767 
0.1752 0.1737 0.1723); 

NAIS@NAIS! +NAIS2: 

ЛћосА15=(Р5А15. /1)%60: Хасе џа! Лос Магі х 

Xplot (NAIS,AoCAIS, '--'), 


J2 


Xthls program computes the takeoff and landing distances for the ANEW aircraft. 
It Is bosed on the anolysis presented In chopter !0 of Hicolol. 


Ulow!85;Xvelocity ot lift off 

T=25400; thrust 

0732. !7; Хоссе! его Гоп дие to grovity 
1253000; Xselght 

Cdo*.02;Xporas!t!c drog 

57639; Xtotal в!пд oreo 
n0*».0023769;Xdens!ty (90 deg. doy--5.00224!) 
C1-2.04:Xcoefflclent of lift 

b*72;Xseing spon 

helf1.4;Xhelght of wlng oboue ground 
Ph*((16*h/b)^2)/(1*((16*h/b)^2)); 
Rn-8.11:Xospect гої їо 

ew. 8: Xeffictency 

K*1/(pl*e*RR); 

L».5*h0*Ulo^2*S*CI;X11 ft 
Cd-Cdo*(Ph*CI^2*K);Xcoefflclent of drog 

02 .5*R0*Ulo^2*5*Cd;Xdrog 

fr*.04:Xfrlction 
$10-(Ulo^2*(1/g))/(2*(T-(D* fr*(H-L)))), Xdlstonce to tokeoff 
Sro*3*Ulo,Xdlstonce to rotote 
Rf*Ulo^2/(g*(1.152-1));Xrod!us of rotat!on 
Sc! -Rf*«1n(.1697B), 
Htof*Rf*(1-cos(.16978)), 
Sobs*(50-11tof)/tan(.16978), 

6tot *S51o*Sro*Scl*Sobs, 
6100*1.414*1^2/(g*R0*S*3*(T-(D* fr*(H-L)))), 


11247000; 

Clm=3; 

Usesaqrt (2*UI/(Cim*AO*S)): 

Ul=1.2*Us; 

Ul f*1.235*Us; 

CIf-2*U41/(RO*UI f^2*5): 

Cd-Cdo*(Ph*CIlm^2*K); Xcoefficlent of drog 
02.5*R0*UI^2*S*Cd;Xdroq 

(г!”-,5: 

RIf=UIf22/(g*(!1.22-1)), 
50!=(50-(В!{*(1-соз(2*#р!/180))))/%ап(2*р!/180), 
SIf=RIfteln(2*pl/180), 
51-1.69*11^2/(g*R0*S*Cla*(T-(D*fr1*(H-L)))), Xlonding rollout 
Stft=Sgl+Stf+St, 
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APPENDIX | 


XThis program sili compute the stabl!itu deríivatiues for thres filght 
conditions. The conditions wil! be at li*0.2, 0.40, 0.76. Corresponding altitudes 
wilt be hes!, 30K, and 30K respectively. These conditlons wll! be denoted by a 
t, 2, and 3 respectiudly. When parameters have defined with little more than an 
educated guess, It wll! be denoted elth a * symbol. Calculatlons are done JAY 
Roskam Part UI. 


M=47100D;Xm!d ronqs welght 
$*639;Xelng reference area 
Lc4*17.5*p1/180: Xsweep at quarter chord 
K*1/(plI*.8*8.11); 
Cdo*0.02;Xparas!t lc drag coefficient 
Caowf=-.1542:%Roskam Part UI, Chap 8 
дСтдс1=-.245; 1 (оСа/оС! )аџегаде ог ОабСот & Roskam results 
0= 1: Хсоипі ег 
Гог П». 2:.28:, 77, 

If «0.3, 

Pe2116.2; Xpressure @ sea level 

else 

P«2116.2*.2975; Xpressure e 30K 

end 

ПП(0)-Н; 

CL(Q)-U*2/(1.4*P*t1^2*5) ; Xcoefflclent of Tift 
Ca(Q)-Cnosf*CL(Q)*dCadCI;X1inear moment coefficient 
CD(Q)-Cdo*K*CL(Q)^2;Xdrag coefflclent 

CDu(Q)-(-4) *K*CL(Q)^2; Xeqn(10.10) 
CLu(Q)-(H^2*coe(Lc4)^2*CL(Q))/(1-H^2*cos(Lc4)^2); Xeqn(10.11) 
0*0*1!:Xcounter 
end 


CLa-[4.822 5.17 6.25]; computed In the LIft Curve Stope progran. 
Caa*dCsadC1.*CLa:Xeqn(10.19) 


Sh*180; Xhonr!zontal! tall surface area 

Kbach*(25.7/9.77);Xdsf!ned !n chapter 10, Page 380 

4bcg*(5.1/9.77);Xdef!ned In chapter 10, Page 380 
ada=.95;%*horlzontal-to-freestream dynamic pressure (qh/q) 
deda0.33;%*downsash gradient at horizontal tal! (page 272) 

CLah*(3.00 3.35 4.43): 8811 ft curve sloped of the horizontal &vertical tolls 
Ubh* (Kbach-Xbcg)*(Sh/S) ;KhorlIzontal tall volume coefficient 
CLad-2*ada*deda*Ubh., *CLah; XC1 alpha dot 

Caad-(-2)*ada*deda*Ubh* (Wbach-Mbcq) . *CLah; Cm alpha dot 


ХТћ! а concludes the longitudinal calculat!ons FOR HOI! and begins Lat-Dir 
Xcalculat lone. 


11) Cuf-sldsforce-due-to-sidssiip (10.2.4.1.1) 
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Dih=2;%dthedral (In degrees) 

Kiet. 73;8from figure 10.0 (2x7-23.5 & df/271) 
Ro=3.5;%radius of fuselage ehere the flow ceases to be a potential (figi0.10, 11) 
9o*pfl*Ro^2:Xarea at that polnt 

Bu*t10:Xtotal span of the vertical tal! 

Su*15;Xarea of one of the vertical talls 

flue*Bu^2/Su: Xvert Ica! tall aspect ratlo 

flurat10*1.028; Xfron figure 10.19 

flueff*Ru*Rurat lo:Xeffect iue Rv 

CyBveff=3;%from figure 10.18 

Суга! 10=0.865;Х from figure 10.17 

CyBer-.00573*DIh;XCyB of the sing 

CyBf=(-2)*KI*(So/S):%CyB of the fuselage 

CyBu=(-2)*Cyrat lo*CyBueff*(Su/S);XCyB of the vertical tall 


12) CIP-rolling nosent-due-to-sldes!tp (10.2.41.1.2) 

CIfCI--.001!;Xfrom figure 10.20. Iterating between taper ratio of OR .5 
KaL=(1.01 1.125 1.3);%ftgure 10.2! using H=.2,.48,.76 & c/2=15 degrees 
Kf-0.97;¥flgure 10.22 

CIBCIA=.0002;Xfigure 10.23 

CIBOIh=-.00022;Xflgure 10.24. Iterating between taper ratio of OR .5 
B-72;Xelng span 

RR-B.11!:Xaspect ratlo 

Ofave=((p!*3.75°2)/. 7854)*.5; 

aC IPO th=(-.0005)*AR* (Dfave/B)*2; 

KmOth=(1.01 1.07 1.2);X%ftgure 10.25 using l1=.2,.48,.76 & c/2*15 degrees 
се=-3.5;Хзев (!диге 10.9 

AC1fze*.012*RR^.5*(2«/B)*(Üfauve/B); 

etan*0.94:X*tan(17.5)tImes wing teist of (-3) degrees. see page 397 
ACIfet--.00003!:Xf1gure 10.26 

for Q=1:3, 
C1fsf(Q)-57.3*(CL(Q)*(CIBCI*KmL(Q) *K f «C TPCIR)*DIh*(CIBDIh*KmDIh(Q)*ACIDDIh) *aCIDz 
e*etan*aCIifiet ); XCIf of the wing-fuselage combinat lon 

end 

Bh*24:Xhortzontal tal! span 

Clfüihf*.65.*CIlfwf; X*CID of the tall-fuselage comblinat lon 
CIBh=(Sh*Bh/(S*B)). *CIBhf;XCIB of the horizontal tall 

2u74;Xsee figure 10.27 

0-24; Хзее figure 10.27 

alfep!/!80*({10 4 0);Xestimated A.0.A from the respective Ci's 
CiPu=CyB*( (Zu. *coalalf)-Lu. *ain(alf))/B);XCIB of the vertical tal! 
CIP-CIPsf+CIBh+CiBu;sXthe grand total 


1X3) Cnfi-yaeIng moment-due-to-sidesiip (10.2.4.1.3) 


Cnfie*0:Xapprox!nmate 
Kne.00165:Xflaure 10.28 
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Kreis (250; 5 tiligurée stone? 

Sfs*376;Xapproximate fuselage side area 

Lf=55;Xfuselage length 

CnBf-(-57.3)*Kn*krI1*(Sfs*L f/(S*B)); XCnB of the fuselage 
CnBv-(-Cyfiv)*((Lu.*cos(al f):2o.*sIn(al f )) /B) ; CnB of the vertical tall 
Cnfi-CnfBe*CnB f *CnBv; Xthe grand total 


X4) CyBd-sideforce-due-to-rate of-sideslip (10.2.5.1) 

Stgba=(-.023 -.025 -.028];Xflqure 10.30 

SIgbd-[.84.87.90];Xflgure 10.31 

Slgbet-[-.02 -.022 -.021];Xflgure 10.32 

SIgbsf*[.14 .145 .15]; XfIgure 10.33 

et*(-3);X*wIng twist in degrees 

Lp*26:Xquarter chord of wing to quarter chord of vertical tall 

2p7l0;Xfrom bottom of fuselage to quarter chord of the vertical tall 
for Oet:3, 

dS1gdfi(Q)-SIgba(Q)*al f(0)*180/pl*Stgbd(Q)*(DIh/57.3)-Slgbet (Q)*et *Slabwf(0);Xean. 
10.47 

CuBd(0)-2*dSlgqdf (Q)*(Su/S) *((Lp*cos(al f(Q)) *2p*sin(al f(Q)))/B);Xean. 10.46 


X5) ClBd-roli!ing monent-due-to-rate of-sidesiip (10.2.5.2) 
Clfid(Q) -CyBd(Q)*((2p*cos(alf(Q))-Lp*sIn(al f(Q)))/B);Xeqn. 10.48 


46) Cnüd-uawing noment-due-to-rate of-sidesitp (10.2.5.3) 
CnfBid(Q) -CuBd(Q)*((Lp*cos(al f(Q))*2p*sIn(alf( Q)))/B);Xean. 10.49 


X7) Cup- sldeforce-due-to-roll rate (10.2.6.1) 
Cup(Q)*2*Cufiv*((Zu*cos(alf(Q))-Lu*sIn(alf(Q)))/B) ;Keqan. 10.50 
end 


X8) Clp- rolling moment-due-to-rol! rate (10.2.6.2) 

for Q*1:3, 

BHa(Q)-*(1-nf(Q)^2)^.5;Xeqn. 10.53 

KHa(Q)*(CLa(Q)*Brta( Q0))/(2*p1);Xean.10.54 

end 

CLaratlo*1;X11ft сое Ис! еп га Го 

BCIpk=s[-.49 -.40 -.43];Xflgure 10.35 
Clpdr=1-4*2Zw/(B*tsin(2*pl/1860))+1I2*(Z/B)“2*(sin(2*pl/180))“2;Xeqn. 10.55 
CipDCLr=-.0015;%flgure 10.36 

COow=.0059;Xfrom the CDo program 

Clph*0;Xapproximate from eqn. 10.59 

Clpu*CyBu*2*(2u0/B)^2;Xeqn 10.60 

Гог 0=1:3, 

Cipdrag(Q)=CipDCLr*CL(Q)*2-.125*CDow;Xeqn. 10.56 
Clps(Q)-BCIpk(Q)*(KHa(Q)/Bita(Q))*CLaratlo*Clpdr*Clpdrag(Q);Xean. 10.52 
end 


Clp*Ciph*Cilpu*Cipwiithe arand total (Тіпе100) 
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X9) Cnp- gyasing mosent-due-to-rol! rate (10.2.5.3) 

Cbar*9,77?; XH.R.C. 

Mbar-0;Xdlstance from the c.g. to the a.c. (posltue for a.c. aft of c.g.) 
Cnpet*.0004:Xf1gure 10.37 
C0*cos(Lc4):C02-*(cos(Lc4))^2; Tfi» tan(Lce4); Tf2-tan(Lc4)^2; 

СпрС100= (1/6) * (АА+6* (АА+ СО) * ( (ХЬаг/СБаг ) * ТАПА +ТА2/12) ) / (АП + 4+0): Хелп. 10.65 
Гог 0=!:3, 

Bnp(0)*(1-111(0)^2*C02)^.5;Xean. 10.64 

CnpCl0lt(Q) *( (RR« 4*C0)/(RR*Bnp(Q)*4*CO)) * (RR*Bnp(Q)*.5*(RR*Bnp(Q)*CO)* TR2) /(RP*.5 
*(RR*C0)*TR2)) *CnpC100;Xeqn. 10.63 

Cnps(Q)-*(-CnpCl0M(Q))*CL(Q)*Cnpet*et:Xean. 10.62 

Cnpu(Q)*(-(2/(B^2))) *Cufiuf(Lu*cos(al f (Q)) *2u*sIn(al f(Q)))*(2u*cos(alf(Q))-Lou*sln( 
alf(Q))-2u);Xeqn. 10.67 

end 

Cnp*Cnps*Cnpu, Xthe grand total 


X9) Сіа- lI!ft-due-to-pltch rate (10.2.7.2) 

Xe-O;Xflgure 10.39 

for Q=1:3, 

Clqwt0(Q)=(.5+2*ke/Cbhar)*CLa(Q):Xeqn. 10.71 

Clae(Q) -((RR*2*C0)/(RR*Bnp(Q)*2*C0))*Claslt0(Q) ; tean. 10.70 
Clah(0)-2*CLah(Q)*Ubh*ada:Xean. 10.72 


end 

Cla*Clqs*Clqh,Xthe grand total 

X И Ее и Е п CO. io ee en mie ms s Mi eem em -- om m ma a 
X10) Caq- pitching moment-due-to-pitch rate (10.2.7.3) 

for Q*1:3, 


Cmq(Q)=1.1*(-2)*CLah(Q)*ada*Ubh*(Xbach-Xbceg) ;¥eqn. 10.78 times 1.1 to account 
Xfor the wlng-body component.This Is from Roskam's "Airplane Flight Dynamics and 
KAutomatic Filght Controls" book Part |, page 188. 


end 

ee ee ctc c uut —AM 
Xback to the lat-der derlvat lues briefly 
ЕВ. жғне Е______.___------------- 
X11) Cyr- sldeforce-due-to-yaw rate (10.2.8.1) 

for (71:3, 


Cur(0)*(-2)*CuyBu*(Lu*cos(alf(Q))*Zu*sln(alf(Q)))/B; ean. 10.80 
end 


%12) Cir- rolling moment -due-to-yaw rate (10.2.8.2) 
CIrCL00-.257;Xflgure 10.41 

ACI rdth=.083*%p 1 "ПА “а | п (6 24)/ (ПА + 400): Хеап. 10,84 
ACIret-(-.014);Xflgure 10.42 

fon omi. 


a 


HUt=1+((RA*(1-Bnp(Q)*2))/(2*Bnp(0)* (AR*Bnp(Q) +2*CO) )) +((AR*Bnp(Q)+2*CO)/(AAR*Bnp(0 
)*4*C0))*TR2/8;Xnumerator of eqn. 10.83 

DE1*1*((RR*^2*CO)/(RR*4*CO))*TR2/B; denominator of eqn. 10.83 
CirCLON(Q)=(HUI/DE!)*CIrCLOO;%eqn. 10.83 
Clre(Q)*CL(Q)*CIrCLOM(Q)*aCInrdih*Dih*aCiret*et:Xean. 10.82 

CIrv(Q)*(-(2/(B^2))) *Cuf'u* (Lu*cos(at f (Q)) *2u*sln(al f(Q))) *(2Zu*cos(alf(Q))-Lu*sin( 
alf(Q)));Xeqn. 10.87 

end 
Cir=Cire+ClrusXthe grand total 


113) Cnr- yawlng monent-due-to-yas rate (10.2.8.3) 

CnrCLr*0;Xflgure 10.44 

CnrCDo*(-.35);Xflgure 10.45 

for Q*1:3, 

Cnre(Q)*CnrCLr*CL(Q)^2*CnrCDo*CDos; Xeqn. 10.87 
Cnrv(0)*(2/(B^2))*Cyfiv* (Lu*cos(al f(Q)) *2u*sIn(al f(Q)))^2;Xean. 10.88 
end 

Cnr-Cnrs*Cnru;Xthe grand total 


Kbe.47;Xflgure 8.52 

CldCidt=.62;%*figure 6.15. Hote:the elevator-to-hor. tall chord ratio & the 
falleron-to-chord rat!o are about thá same. Thís 19 Important for section 17). 
C1dt*5.2; figure 8.14 

Kprise*!;Xapproximate (figure 8.13) 

RdCLRdc!-1.02; figure 8.53 
RIfde-Kb*CldCidt*Cidt*RdCLRdci*(Kprine/(2*p1*.88));X*eqn. 10.94 


1314) Clae- lift-due-to-elevator (10.3.2.2) 
Гог 0=1:3, 
CL!h(Q)-ada*(Sh/S)*CLah(Q);Xeqn. 10.91 
Clae(Q)-Ri fde*CLIh(Q):Xean. 10.95 

end 


X15) Cnae- pitching aoment-due-to-eleuator (10.3.2.3) 
for (91:3, 

Calh(Q)*ada*Ubh*(-CLah(Q)):Xeqn. 10.91 

Cmae(Q)*Al fdatCaih(Q):eqn. 10.95 

end 


116) Cyaa- sideforce-due-to-alleron (10.3.5.1) 
Cy4a*0;Xeqn. 10.105 


117) Claa- rollIna nosaent-due-to-alléron (10.3.5.1) 
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bCplak-[.4 .395 .385];Xflqure 10.46b 

ог ОЗЕ "3. 

Cpla(Q)-(Kfa(Q)/Blla(Q)) *bCptak(Q);Xean. 10.107 
Alfdela(Q)=(Cidlidt*Cidt)/CLa(Q);Xeqn. 10.109. 
Cla(Q)-Rifdela(Q)*Cpla(Q):Xean. 10.108 

end 

Claae2*Cla:Xeqn. 10.113 


118) Cnaa- yawing moment-due-to-alleron (10.3.5.1) 
Kas-.115;Xflgure 10.418 

(ог 0=1:3, 

Cnaa(Q)*Ka*CL(Q)*Claa(Q):Xean. 10.114 

end 


1319) Cuar- slideforce-due-to-alleron (10.3.8.1) 
$u2790:Xtotal! vertical tall area 

Kp2*.8;Xflgure 8.13 

CldCldt2-.82;X*flgure 8.15 

Cldt2-5.7;Xflgure 8.14 

for Q*1:3, 
Cyar(Q)-CLah(Q)*Kp2*Kb*CldCidt2*Cidt2*(Su2/S) ; Xeqn. 10.123 
end 


20) Clar- rolling sosent-due-to-alleron (10.3.8.2) 

for Q*1:3, 

Clar(Q)-Cuyar(Q)*((Zu*cos(al f(Q))-Lu*sIn(alf(Q)))/B);Xean. 10.124 
end 


X21) Cnar- gawing mosent-due-to-alleron (10.3,8.3) 

for Q*1:3, 
Cnar(Q)*(-Cyar(Q))*((Lu*cos(atf(Q))*2u*s1n(at1f(0)))/8):Xeqn. 10.125 
end 
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APPENDIX J 


XThis program sil! calculate the dynomic characteristics of the AEH aircraft. 
The programming !s based on the dynamic approximations presented In Etkin's 
baak, First edition, 31959, Chapters 6& 7. Stability Derivatives are acquired 
from the Stablilty бег уа ісе progran. 


Hass*53000/32.2:%mass In sluga 

Сбаг=9, 77; Хаеап аегадупат!с chord 

5639; Хе! пад reference area 

Li*Cbar/2;Xpage !92 (langltud!nal only) 

RD!*.D0023769; Xdensltu at sea level 

R02*.0023769*.31D6; Xdensltu at 350DD ft. 

НИ! =Наза/(101#5% | ); Храде 192 

IlU2-fass/(R02*S*L1); Xpage 192 

CL*[1.24113 0.7244 0.2890]: Xreference CL. From Stab. Der. program 
С0=[0.0956 0. 0457 0.024! ]; Xreférence CD. Fram Stab. Der. program 
CLa-[4.8220 5.1700 6.2500]; reference CLa. Frame Stab. бег, ргадгат 
CDu-[-D.3024  -0.1030  -0.0164];Xreference CDu. Fram Stab. Der. program 
alf*pl/180*[10 4 0]:XeetIeated R.D.fi from the reepectiue Cl's 
yc c RN ес а а КК ee COO EE 
Xphugold modes 

Hnp(1)=CU(1)/(sqret(2)%MU1);Xean.(6.7,4) assuming negligible Czu and Cza 
inp (2)"CL(2)/(sart (2) *NU2):Xeqn.(6.7,4) assuming negligible Czu and Cza 
Wnp(3)=CL(3)/ (sqrt (2)*NU2): Хеап. (6.7,4) assuming negligible Czu and Cza 
far Q=1:3, 

Cxu(0)*(-2)*(CD(Q)*CL(Q)*tan(alf(Q)))-COu(Q) ;Xpage 150 (11) 

Zep(Q) *(-Cxu(0))/(2*sqrt (2)*CL(Q)) ;eqn.(6.7,4) assuming negligible Czu and Czq 
ldp(Q)*sart(1-2ep(0)^2)*lInp(Q) ; Xdaaping frequency 
To(Q)=(2*pl)/Ndp(Q) ; Xperlod 

end 

Charíi*[1 (2*2ep(t1)*tnp(1)) Wnp(1)^2];Xcharacter!st!c equat lon 

Сћаг2= [1 (2*2ep(2)*Unp(2)) !Inp(2)^2];Xcharacterlstlc equat lon 

Char3-[! (2*2ep(3)*tnp(3)) Unp(3)^2];Xcharacterlst Ic equat lon 
Al=roats(Char!);%the raots 

R2-raote(Char2):Xthe raats 

R3*raats(Char3):;Xthe roots 


Xshort perlad mades 

10у 74176; Хтатепі аї (пегі 1а {гат the CG pragran 

Ibi*Iyg/(R01*S*L1^3); Xnon-dImenslona! mament of Inertla. Page 192, 
Ib2-1yy/(R02*S*L1^3) ; Xnon-dImenslonal mament af Inertla. Page 192. 
Cza-(-1)*(CLa*CD);Xeqn.(5.2,3) 

Cma-[-1.1814  -1.2666  -1.5312];Xfras stablllty derluvatiue pragram 
Caq-[-7.8521 -8.76082 -11.5949]; from etabllity derivative program 
Свад=(-2.3556 -2.6304 -3.4785);Xfrom stablilty derlvative program 
Uns(1)esart ((Cza(1)*ma( 1) -2*NU1*Cma(1))/(2*NUItIbt)) sXean.(6.7,7) assuming 
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negilgibie Czodot and Cza 

for Qe2:3; 

lins(Q)esqrt ((Cza(Q)*Cmq(Q)-2*NU2*Cma(Q))/(2*NU2*1b2));%eqn.(6.7,7) assuming 
negligtble Czadot ond Cza 

end 

Zes(1)-(-1)*((2*ttUI *Cmna(1)*Ibf*Cza( 1)*2*fUI *Cmad(1))/(2*(2*HU1I* Ib1*(Cza( 1) *Cma(1) 
-2*HnU1*Cna(1)))^.5)):Xe6qn.(6.7,7) assusing nealigibie Czadot and Ста 

for Q*2:3, | 
Zes(Q)«(-1)*((2*HU2*Cma(Q)*lb2*Cza(Q)*2*HtU2*Cmad(Q))/(2*(2*IU2* Ib2*(Cza( 0) *Cma(Q) 
-2*NU2*Cma(Q)))*.5));Xeqn.(6.7,7) assuming negligible Czadot bnd Cza 

end 

for Q*1:3, 

Hde(Q)*sqrt(1-2es(Q)^2)*tIns(Q); Xdanplng frequency 
Ts(Q)*(2*pl)/HUds (Qo) ; Xper lod 

end 

Chartis-[I (2*2es(1)*1ns(1)) Hns(1)^2]:Xchoracteristic equation 

Char2s*[l! (2*269$(2)*1Ins(2)) !Ins(2)^2]; Xcharacterlstlc eauct lon 

Char3s-[(1 (2*2e5$(3)*1Ins(3)) Uns(3)^2]:Xcharacteristic equat lon 
Als=roots(Charls):X{the roots 

R2s"roots(Chor2s) :Xthe roots 

h3se*roots(Char3s) :Xthe roots 


B=72;Xelng span 

L2-*B/2;Xpage 226 

Ixx*100006; X»oment of Inertla from the CG program 

І22"147693:Хтотепі о” Inertia from the CG program 

Ixz=-14.9335; Xmoment of Inertla from the CG program 
lal*Ixx/(R01*5*L2^3) ; Xnon-dImenslonal moment of Inertla. Page 192. 
la2*Ixx/(R02*5*L2^3) :Xnon-dImenslonal moment of Inertla. Page 192. 
Ici*Izz/(h01*9*12^3) ; Xnon-dImenslonal moment of Inertla. Page 192. 
Ic2-I1zz/(R02*5*L 2^3) :Xnon-dImenslonal moment of Inertlo, Раде 192. 
lel2Ixz/(R01*S*L2^3) ; Xnon-dImenslonal monent of Inertia. Page 192. 
le221xz/(h02*5*L 2^3); Xnon-dlImenslonal moment of Inertla. Page 192. 
CyB--0.5877;Xfrom stabllity derivative program 

Cyr=0.2437;Xfrom stablilty derluative program 

Clp-[-2.4765 -2.5993 -2.B8140);Xfrom stability dertvative program 
Сіге(0.4717 0.3620 0.2667);Xfrom stabllity derlvative program 
Cnp*[0.1319 0.0764 0.0291);Xfrom stabllity derluative program 
Спг=[-0.0655 -0.0848 -0.08633);X%from stablilty derluatlve program 
С1В=[-0.1279 -0.1307 -0.1273);%from stabllity derluative program 
Cyp-[0.0023 -0.0235 -0.0406);Xfrom stability derlvative program 
CnB- [0.0576 0.0571 0.0560);%from stobllltu derluatlue program 


R(1)-2*1U1*(laf*Icl-let^2); Xpolynomlal coefficient. ean.(7.1,3) 
R(2)-2*t1U2*(1a2*1c2-1e2^2) :Xpolunomlal coefflclent. ean.(7.1.3) 
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n(3)-R(2); 
B(1)-CuyB*(1el^2-lal*Icfl)-2*HUt*(Ict*Clp(1)*lal*Cnr(I)*lel*(Clnr(1)*Cnp(1)));Xpolun 
omlal coefficient. edan.(7.1,3) 

for (22:3, 

B(Q)-Cufi*(1e2^2-1a2*1c2)-2*tU2*(Ic2*CIp(0)*la2*Cnr(Q)*le2*(Clr(Q) *Cnp(Q)));Xpolun 
omlal coefficient. eqn.(7.1,3) 

end 

C(1)72*tUt * (Cnr (1)*CIp( 1) -Cnp(CI)*CiIr(1)*lal*CnBCT)*1et*CIB(1))*lot *(CuB *Cnr( 1)-Pn 
BCT )*Cyr)+tol* (CyB CI p(l)-CIB(1) *Cyp lt) )+lel *( CyB Cnp (1 )-CnB(1) *Cyp(l )+Cir (1) * CyB 
-Cyr*CIB(1)):Xpolynomlal coefflctent. eqn.(7.1,3) | 

Гог 0“2:3, 
С(0)=2*102*(©пг(0)*С!р(0)-Спр(0)*С!г(0)+!а2*Сп3(0)+1е2*С10(0))+1ва2*(СуВ*Спг (0)-Сп 
0(0)*Cur)+|c2*(CuD*Clp(0)-C1I0(0)*Cun(0))+le2*(Cu n *Cnn(Q)-Cnp(Q)*Cun(0)+CIr(0)*Cuñ 
-C[ur*C10(0));Xpolunonlal coefflclent, eqn.(7.1,3) 

end 
0(1)”СуПВ*(С1г(1)*Спр(!1)-Єпг(1)*С!р(1))+Сур(!)*(С!П(1)*Спг(!)-СлП(1)*С1г(1))+(2*!!! 
I-Cyr)*(CI8(1)*Cnp(T1)-CnB(C T) *CÓp( 12) -CLCT)*CIe 1*CT0(1)* Tet *CnB(1)) ; Xpolunomlal 
coefficient. eqn.(7.1,3) 

Гог 0-2:3, 

0(0) «Суб *(С1г(0) *Спр(0) Сп" (0) *С1р(0) ) +Сур(0)* (С18(0) *Спг(0) -СпВ (0) *С1г(0) ) + (2 * 11 
2-Суг) %(С16(0) *Спр(0) -СпВ (0) *С1р(0) ) -Сі (0) * ( 1с2*С18 (0) + 1е2 Сп (0) ); Хро і упот! а! 
coefficient. вап. (7. 1,3) 

end 

E(1)=CLU(1)*(C1I0(1)*Cnr(1)-CnD(1)*CIr(1));Xpolunomlol coeffictent. eqn.(7.1,3) 

(ог 0=2:3, 

E(Q)-CL(Q)*(CIB(Q)*Cnr(Q)-CnRB(Q)*Clr(Q));Xpoluynomlal coefflctent. eqn.(7.1,3) 

end 


ChoarLD1*[R(1) B(1) CC1) 0(1) E(t)]:Xcharacter!stic equatlon 
ChorL02-[R(2) B(2) C(2) D(2) E(2)]:Xcharacterlst !c equat fon 
CharLD3-*[R(3) B(3) C(3) D(3) E(3))]:Xcharacteristic equat lon 
nLDI*roots(ChorLD1); Xthe roots 

RLD2-*roots(CharLD2) :Xthe roots 

RLD3-roots(CharLD3) : Xthe roots 

[UnL1,2eL1] = DANP(CharLO!);Xnatural frequency and damping rat lo 
[InL2,2eL2] - DfilfP(CharLD2) ;Xnatura! frequency and damping roto 
[unL3,2eL3] * DRIIP(ChanLD3); Xnatural frequency and dompling rot lo 
HdL1-7sqrt(1-2eL1.^2).*tInL 1; Xdamping frequency 
TL1*(2*p1)/HdL 1 ; Xper lod 

HdL2-sqrt(1-28L2.^2).*lInL2; Xdamping frequency 
TL27(2*pl)/HdL2; Xper lod 

HdL3-sqrt (1-2eL3.^2). *IInL3; Xdamplng frequency 
TL3*(2*p1)/lldL 3; Xper lod 
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